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INTRODUCTION 


The  general  purpose  of  this  project  is  to  demonstrate  the  possibility 
of  controlling  some  processes  in  the  ionosphere  by  application  of  radio 
energy  at  the  gy rof requency .  Another  purpose  is  to  estimate  the  magnitude 
of  the  effects  produced,  to  relate  them  to  the  properties  of  the  ionosphere, 
and  thus  eventually  to  devise  new  methods  of  probing  the  ionosphere. 

Another  prupose  of  the  contract  is  to  perform  laboratory  experiments 
on  the  interaction  phenomena  in  an  ionized  gas  subjected  to  a  steady 
magnetic  field,  and  to  scale  the  results  of  these  experiments  to  the  actual 
ionosphere  in  view  of  its  possible  control. 

An  experiment  has  been  planned  and  executed  which  consisted  of  firing 
a  rocket  through  the  D  and  lower  E  regions  of  the  ionosphere.  The  rocket 
carried  a  transmitter  which  periodically  radiated  a  strong  pulse  at  the 
gyrof requency  (about  1.4  me)  The  temperature  of  the  electrons  in  the 
vicinity  of  the  rocket  having  been  raised  by  this  ’heating"  signal,  the 
properties  of  the  medium  (collision  frequency,  rates  of  various  processes, 
and  the  electron  density,  if  the  signal  were  strong  enough)  would  have  been 
modified  These  modifications  were  detected  by  observing  amplitude  modula¬ 
tion  upon  a  sensing  signal  propagated  through  the  heated  region. 

The  planning  phase  of  the  project  began  in  mid  1959  The  analogy 

between  the  rocket-borne  experiment  and  the  Luxembourg  effect'  was  helpful; 

however,  the  differences  were  important  and  necessitated  careful  study  in 

several  areas.  Graphical  methods  were  developed  for  the  computation  of 

12  3 

refractive  index  and  polarization  in  a  magneto  ionic  medium  ’  '  .  Maxwell' s 

equations  were  solved  for  the  case  of  an  infinitesimal  dipole  antenna  in  an 

4 

anisotropic  plasma  and  both  radiation  patterns  and  near  fields  were  obtained  . 

Other  studies  more  directly  related  to  the  rocket  borne  experiment 
have  been  included  in  progress  reports  These  include  the  choice  of  sensing 
wave  frequency,  heating  antenna  system,  and  heating  signal  parameters  (such 
as  pulse  width,  pulse  repetition  frequency  and  power  level) 

Chapter  1  of  this  report  concerns  the  planning,  instrumentation  and  pre¬ 
liminary  analysis  ol  results  obtained  from  the  gy roresonance  probe  carried 
by  Aerobee  sounding  rocket  AA1  193  on  1  May  1962  In  addition  to  the  detection 


2 


I 
I 

of  cross  modulation  on  the  sensing  wave  the  experiment  included  the  measurement  1 

of  the  polarization  ellipse  of  the  sensing  wave,  the  gyrof requency  radiated 
power,  radiated  field  intensity,  and  antenna  impedance.  Nonlinear  effects  « 

measured  included  the  excitation  light  produced  near  the  heating  antennas  I 

and  the  difference  frequency  generated  by  mixing  of  the  heating  and  sensing 
signals  in  the  ionosphere.  System  performance  was  monitored  along  with  the  J 

aspect  of  the  earth  s  magnetic  field  Unfortunately  much  of  this  data  was 

I 

not  obtained  due  to  the  failure  of  motor  drivenswi tch.  This  device  served  J 

the  dual  purpose  of  data  sampling  commutator  and  clock,  its  failure  resulted 
in  loss  of  both  synchronization  and  time  shading  modes  of  data  acquisition. 

However  the  sensing  wave  modulation  and  polarization  ellipse  as  well  as 
gyrof requency  field  intensity  were  obtained  intact  The  analysis  of  this 
data  is  far  from  complete  and  will  be  included  in  work  done  on  the  next 
contract.  This  experiment  has  provided  useful  information  to  guide  the 
instrumentation  of  future  gyroresonance  probes.  It  has  also  demonstrated 
the  applicability  of  the  electron  heating  technique  in  studying  the  composi¬ 
tion  of  the  lower  D  region  by  modulation  of  Faraday  rotation. 

In  conjunction  with  the  planning  of  the  rocket  experiment  it  was 
necessary  to  carry  out  plasma  experiments  in  the  laboratory  and  to  relate 
the  results  to  ionospheric  problems  Included  were  studies  of  the  ion  sheath, 
plasma  heating,  and  antenna  impedance  A  theory  for  the  impedance  of  a 
short  antenna  in  an  anisotropic  plasma  was  developed,  verified  experimentally, 
and  used  to  design  the  heating  antenna  matching  network  for  the  rocket  experi¬ 
ment.  This  work  is  discussed  in  Chapter  2 

The  results  of  a  laboratory  study  of  the  nonlinear  interaction  of  an 

5 

RF  electromagnetic  wave  and  a  plas.na  have  been  published  and  are  summarized 
in  Chapter  3.  The  l nst rumentat l o.i  of  the  difference  frequency  measurement 
for  the  rocket  experiment  vas  included  in  (he  course  of  this  work. 
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1.  GYRORESONANCE  PROBE  EXPERIMENT  -  R.  R.  Hodges 

1.1  Introduction 

Cross-modulation  occurring  in  the  ionosphere  was  first  reported  by 
Tellegen  in  1933.  This  effect  is  commonly  known  as  the  "Luxembourg  effect" 
because  it  was  observed  that  modulation  of  the  Radio  Luxembourg  station 
(250  kc)  was  impressed  upon  weaker  signals  traversing  the  region  of  the 
ionosphere  near  this  station. 

The  first  theoretical  explanation  of  the  "Luxembourg  effect"  was  pre- 

7 

sented  in  1934  by  Bailey  and  Martyn  .  They  evaluated  the  average  energy 
received  by  electrons  between  collisions,  deduced  the  increase  in  collision 
frequency,  and  obtained  a  formula  for  the  cross  modulation.  Refinements 
upon  this  theory  have  been  made  and  much  experimental  work  has  been  done. 

g 

Notably  a  prediction  by  Bailey  in  1938  that  a  resonance  effect  should 

occur  when  the  disturbance  is  produced  by  a  signal  at  the  gyrof requency 

9  10  8 

has  been  verified  ’  An  additional  prediction  by  Bailey  was  that  a 

small  amount  of  power  (1  watt  for  example)  at  gyroresonanee  and  concentrated 

in  a  sufficiently  narrow  beam  should  produce  an  observable  amount  of  cross 

modulation  in  the  lonospherS. 

Due  to  the  impractical  it les  of  the  required  antenna  system  (one  wave¬ 
length  at  the  gyrofrequency  being  about  200  meters),  Bailey  s  latter  pre¬ 
diction  has  not  been  verified  experimentally. 

However,  by  carrying  the  gy roresonant  transmitter  into  the  ionosphere 
with  a  rocket  the  antenna  problem  is  alleviated.  All  of  the  radiated 
gyrofrequency  energy  is  available  to  heat  electrons. 

A  probe  designed  to  perform  such  an  experiment  was  carried  to  an  altitude 
of  94  km  by  Aerobee  sounding  rocket  AA1  193  on  1  May  1962  The  disturbance 
of  the  ionosphere  was  produced  by  a  rocket  borne  gyrofrequency  transmitter. 

It  was  detected  by  a  sensing  wave  radiated  from  the  ground  through  the 
ionosphere  and  received  at  the  rocket  The  cross  modulation  was  telemetered 
and  recorded 

In  addition  to  the  measurement  of  cross  modulation  the  probe  was 
instrumented  to  obtain  data  from  which  the  gyroresonanee  heating  antenna 
impedance  and  radiated  power  could  be  obtained.  The  near  electric  field 
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produced  by  the  heating  signal  and  the  excitation  light  due  to  inelastic 
electron-molecule  (dr  atom)  collisions  near  the  heating  antennas  were  also 
measured „ 

The  measurement  of  nonlinearities  in  the  ionospheric  properties  was  pro¬ 
vided  by  a  receiver  tuned  to  the  difference  frequency  of  the  sensing  and 
heating  signals. 

Monitoring  of  the  operation  of  the  system  was  accomplished  by  means 
of  a  data  sampling  switch  Such  things  as  magnetic  aspect  and  sensing 
receiver  AGC  voltage  were  included 

The  data  sampling  function  occupied  one  deck  of  a  two  pole  switch. 

The  other  was  used  as  a  clock  to  provide  synchronization  of  the  experiment. 
This  device  failed  mechanically  at  burnout.  As  a  result  the  data  was  limited 
to  the  cross  modulation  experiment,  the  heating  field  intensity,  and  the 
sensing  receiver  AGC  voltage.  (The  latter  is  the  point  at  which  the  data 
sampling  switch  fortuitously  happened  to  stop.) 

The  cross  modulation  data  does  exhibit  the  expected  results.  However 
the  lack  of  subsidiary  measurements  limits  the  accuracy  of  interpretation 
of  this  data. 

In  subsequent  sections  of  this  part  of  the  report  the  planning  of  a 
gyroresonance  probe  experiment,  the  instrumentation  of  an  initial  experiment, 
and  a  preliminary  analysis  of  the  results  of  this  experiment  are  given. 

1.2  Properties  of _ the  Lower  E  Region  of  the  Ionosphere 

The  choice  of  experimental  parameters  for  the  gyroresonance  probe  was 
based  upon  a  set  of  assumptions  concerning  the  properties  of  the  lower  E 
region  of  the  ionosphere.  These  properties  are  as  follows 

(1)  The  E  region  is  a  weakly  ionized  gas  which  means  that  the  electron 
number  density,  ne,  is  much  smaller  than  the  number  density  of 
molecules,  N  For  example 

at  70  km  n  =  1  5  \  108  m~3  while  N  =■  2  1  X  1021  m'3 

e  11-3  lft't 

at  100  km  n&  =  10  m  while  N-78\10m 

(Electron  densities  are  from  Way n l ck^  and  molecular  densities 
1 2 

from  Minzner,  et  al .  ) 

(2)  The  gas  is  electrically  neutral  on  the  average,  and  negative  ions 
are  assumed  to  have  an  inconsequential  density  Some  fluctuations 
of  net  charge  may  occur  over  regions  of  the  order  of  the  Debye 
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length  X^,  This  length,  given  by  6 . 9 z/x/n  cm,  can  be  a  few  centi¬ 
meters  in  the  lower  E  region  and  becomes  much  less  at  higher 

o  12 

altitudes.  At  100  km  the  temperature  is  200  K  ,  the  electron 

4 

thermal  velocity  is  9.5  X  10  m/sec  and  which  can  be  inter¬ 
preted  as  the  distance  traveled  with  this  velocity  during  a  radian 

7 

period  of  the  plasma  frequency  (t~  ^  2  X  10  ),  is  only  4.8  mm. 

8—3  5 

At  70  km,  taking  n  =  1.5  x  10  m  ,  hence  u/  =  7  X  10  and  T  = 
e  ,  ’  p 

o  4 

210  K,  hence  v  =  9  8  X  10  m/sec,  the  Debye  length  is  about 
10  cm.  Thus  the  Debye  length  is  small  compared  to  the  dimensions 
of  the  rocket  or  the  wavelength  for  the  gy rof requency  (approximately 
200  meters). 

(3)  The  collisions  of  electrons  occur  mainly  with  neutral  molecules. 

This  is  a  consequence  of  (1)  and  (2). 

(4)  For  low  enough  energies  electron-molecule  collisions  are  elastic. 

Since  the  molecules  are  much  heavier  than  the  electrons,  a  collision 
produces  only  a  small  energy  transfer  and  mostly  a  change  in 
momentum.  For  a  perfectly  elastic  collision  between  an  electron 

of  mass  m  and  a  molecule  of  mass  M  the  fraction  of  energy  lost  per 
collision  is  given  by 

G  -  't  2m 

m  *  M  M 

The  mean  molecular  weight  of  air  is  28.766  at  sea  leveal  and  28.90 

12  c 

at  100  km  ,  therefore  taking  M  Q  29,  G  =  4  X  10  ,  When  the  col¬ 

lisions  are  not  elastic  (for  higher  average  electron  energies) 
we  shall  still  assume  that  the  fractional  loss  of  energy  per  collision 

_  3 

is  represented  by  a  constant  G  =  10 

(5)  The  molecules  move  very  slowly  compared  to  the  electrons  At  thermal 

equilibrium  the  ratio  of  the  rms  velocities  iSz/M/m  Q  250  [At 

100  km,  if  T  ~  200  K  v  ^  9.5  X  10^  meters  per  second  v  =?  380 

e  '  m 

meters  per  second  ]  The  velocity  of  molecules  may  be  neglected 
in  the  sense  that  their  distribution  and  the  manner  in  which  it 
changes  with  electron  temperature  are  not  important  They  form 
a  background  capable  of  taking  energy  from  the  electrons  at  a 
certain  rate  and  capable  of  providing  changes  of  electron  momentum 
at  each  collision 


I 

I 

(6)  When  an  electromagnetic  1'ield  is  applied  to  the  medium  it  will  act  I 

mostly  on  the  electrons,  changing  their  momentum  The  collisions  * 

have  two  effects: 

(a)  They  " thermal ize'  the  kinetic  energy  of  the  electrons,  i.e. 
orient  the  velocities  randomly 

(b)  By  transferring  a  fraction  of  the  energy  to  the  molecules  they 
prevent  the  electron  energy  from  increasing  indefinitely.  Since 
molecules  are  much  more  numerous  than  electrons  they  are  capable  of 
taking  this  energy  without  appreciably  changing  their  statistical 
properties  (for  example  the  molecular  temperature  or  cross  section 
for  electron  collision)  To  be  sure  some  molecules  will  be  excited 
or  even  ionized  by  non-elastic  collisions  but  the  proportion  of 
those  that  are  affected  will  remain  small.  The  assumption  of 

G  >  2m/M  in  (4)  above  takes  this  into  account 

(7)  The  scattering  of  the  electrons  by  the  molecules  may  have  some 
directive  properties.  Lacking  information  on  this  point  we  shall 
assume  the  scattering  to  be  isotropic.  This  means  it  is  unnecessary 
to  distinguish  between  the  frequency  of  collisions  and  the  collision 
frequency  for  transfer  of  momentum. 

(8)  The  effective  cross  section  can  be  expected  to  depend  on  the 
electron  velocity  relative  to  the  scatterer  In  the  first  approxi¬ 
mation  it  is  sometimes  assumed  to  be  constant  Then  the  collision 
frequency  v  will  be  proportional  to  the  square  root  of  electron 
temperature  and  to  the  number  density  of  molecules.  When  the  electron 
temperature  differs  from  the  molecular  temperature  we  shall  assume 


o  o 


where  is  the  value  of  collision  frequency  in  thermal  equilibrium 
fi  e  T^  r  T^) ,  and  v  and  v  are  rms  electron  velocities  for  heated 
and  equilibrium  cases  respectively 

(9)  Relevant  parameters  that  are  functions  of  altitude  are  temperature, 
electron  density  and  collision  frequency  These  are  plotted  in 
Figure  1.1  according  to  the  following  sources- 
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(a)  Temperature  from  Minzner,  et  al . 
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(b)  Electron  density  (noon  average)  from  Waymck  and  collision 
frequency  from  Nicolet^V 

1 . 3  Experimental  Procedure 

The  experiment  designed  to  study  electron  gy roresonance  heating  within 
the  ionosphere  uses  a  rocket-borne  pulsed  gyrofrequency  transmitter  as  the 
source  of  disturbance  of  the  medium  The  detection  of  electron  heating  is 
performed  in  the  same  manner  as  that  commonly  used  in  investigating  the 
"Luxembourg  effect",  i.e.  a  continuous  sensing  wave  is  propagated  through  the 
heated  region.  Cross  modulation  detected  on  the  sensing  wave  indicates  the 
extent  of  the  heated  region  weighted  by  the  degree  of  heating  along  the  path 
of  the  sensing  wave. 

Figure  1.2  shows  a  schematic  of  the  cross  modulation  experiment.  The 
sensing  wave  propagates  from  a  ground  based  transmitter  to  a  receiver  in  the 
rocket;  this  insures  that  it  is  sufficiently  weak  so  as  to  have  negligible 
effect  upon  electron  energy.  The  direction  of  propagation  of  the  sensing 
wave  and  the  geomagnetic  field  are  near  enough  to  being  parallel  to  justify 
use  of  the  quasi-longitudinal  approximations  for  the  Appel ton-Hartree  equations, 
provided  the  experiment  is  performed  at  Eglin  Air  Force  Base,  Florida.  Note 
also  that  the  portion  of  the  heated  region  through  which  the  sensing  wave 
passes  is  that  heated  by  a  quasi-longitudinal  gyrofrequency  wave  in  the 
first  approximation  The  output  of  the  sensing  receiver  is  telemetered 
using  a  wide  bandwidth  channel  and  recorded  on  magnetic  tape  at  the  APGC 
telemetry  facility 

The  size  of  the  heated  region  is  limited  by  the  absorption  of  the  gyro¬ 
frequency  wave  as  it  transfers  energy  to  the  electrons.  For  example  the 
attenuation  of  a  gyrofrequency  wave  of  extraordinary  sense  of  polarization 
propagating  parallel  with  the  magnetic  field  is  approximately  3  db  per 
meter  at  85  km  and  8  db  per  meter  at  92  km  prior  to  heating  The  effect  of 
electron  heating  is  to  decrease  absorption  of  this  wave  which  results  in  a 
gradual  increase  of  the  disturbed  region,  but  not  to  dimensions  greater  than 
the  order  of  one  kilometer  This  subject  is  discussed  further  in  Section  1.5. 

As  a  result  of  the  limited  region  of  gy roresonance  heating  the  experiment 
may  be  considered  to  occur  in  a  nearly  homogeneous  medium  The  local  values 
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Figure  1.2.  Gyroresonance  probe-  experimental  schematic 
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of  electron  density,  collision  frequency  and  molecular  temperature,  which 
are  functions  of  altitude,  are  variable  parameters  in  the  experiment.  In 
addition  the  peak  power  output  of  the  heating  transmitter  is  continuously 
varied  through  three  power  levels  separated  by  10  db  to  provide  another 
parameter. 

The  altitude  range  of  the  experiment  is  determined  primarily  by  the 
gy roresonance  bandwidth,  which  is  deiined  as  v/271*  Below  76  km  there  is 
no  resonance.  At  88  km  the  fractional  bandwidth  is  10$>  At  100  km  it  is 
\°jo ,  However,  the  accuracy  with  which  the  gyrof  requency  as  a  function  of 
altitude  can  be  predicted  is  roughly  -  1%  for  a  dipole  approximation  of 
the  geomagnetic  field.  Therefore  the  rocket  flight  should  be  limited  to  a 
maximum  of  100  km  to  insure  that  the  heating  signal  frequency  is  within 
the  resonance  and  that  the  probe  remains  in  the  desired  part  of  the  iono¬ 
sphere  as  long  as  possible. 

The  geomagnetic  field  at  Eglin  Air  Force  Base  is  0.5276  gauss.  Based 
on  the  inverse  cube  law  the  gy rofrequency  at  100  km  is  1.41  me  This  fre¬ 
quency  is  within  the  resonance  bandwidth  for  all  lower  altitudes. 

In  the  instrumentation  of  the  initial  gyroresonance  probe  several  sub¬ 
ordinate  experiments  were  performed.  Included  with  the  heating  transmitter 
subsystem  were  the  instantaneous  measurement  of  heating  antenna  impedance, 
power  radiated,  and  radiated  near  electric  field.  These  were  intended  as 
indicators  of  design  improvements  for  future  experiments  as  well  as  data 
relevant  to  interpretation  of  simultaneous  cross  modulation  measurements. 

A  10  stage  photomultiplier  (type  7767)  was  included  to  detect  excitation 
light  induced  by  inelastic  collisions  of  heated  electrons  and  molecules 
(or  atoms) 

In  conjunction  with  the  investigation  of  non-linearities  in  plasmas  by 
J.  Verdeyen  (see  Chapter  3)  a  receiver  to  detect  mixing  of  the  sensing 
(2,02  me)  and  heating  (1  41  me)  signals  was  included  The  receiver  was 
tuned  to  the  difference  frequency  (610  kc), 

A  functional  block  diagram  of  'lie  entire  gyroresonance  probe  experiment 
is  shown  in  Figure  1  3 

The  block  denoted  System  Operation  Monitors  includes  the  sampling  of 
voltages  which  indicate  the  status  of  various  components  (such  as  deployment 
of  whip  antennas  power  -.upply  operation  sensing  receiver  AGC,  and  magnetic 
aspect ) 
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Telemetry,  data  conditioning,  and  synchronization  will  be  discussed 
later . 

1.4  Probe  Configuration 

The  gyroresonance  probe  instrumentation  was  installed  in  a  fiberglass- 
epoxy  ogive  nose  cone  and  a  15  inch  extension  section.  Telemetry  and  range 
safety  functions  utilized  two  10  inch  extensions. 

Figure  1.4  shows  the  locations  and  orientations  of  the  antennas 
associated  with  the  experiment.  Telemetry  and  beacon  antennas  are  omitted. 

The  gyroresonance  heating  antennas  consisted  of  two  base  loaded 
fiberglass  whips,  each  10  feet  long.  During  launch  they  were  folded  down 
along  the  rocket  body  and  were  to  be  deployed  at  approximately  70  km. 

Loop  antennas  located  within  the  ogive  served  the  sensing  wave  and  dif¬ 
ference  frequency  receivers.  The  sensing  wave  loop  antenna  and  the  heating 
antennas  were  orthogonally  oriented  to  reduce  coupling  of  the  heating  signal 
into  the  sensing  receiver. 

The  electric  field  intensity  of  the  heating  signal  was  detected  by 
measuring  the  RF  voltage  between  two  parallel  plates  located  within  the  nose 
cone  near  its  upper  extremity  (approximately  2  meters  from  the  heating 
antennas) . 

Excitation  light  was  detected  by  means  of  a  photomultiplier  tube  at 
one  end  of  a  10  inch  long  collimating  pipe.  The  other  end  terminated  in  a 
window  at  the  skin  of  the  rocket.  Its  axis  was  oriented  to  intersect  one  of 
the  deployed  heating  antennas  at  a  point  8  feet  from  the  rocket.  The  3  db 
beamwidth  was  approximately  3°. 

The  coordinates  shown  in  Figure  1.4  are  used  to  designate  the  axes  of 
three  magnetic  aspect  sensors.  The  use  of  three  magnetometers  gives  directly 
the  magnetic  aspect  of  heating  (x)  and  sensing  (y)  antennas  as  well  as 
precession  of  the  z  axis 

1.5  Gyroresonant  Heating  of  Electrons  within  the  Ionosphere 

A  rigorous  treatment  of  the  problem  of  electron  gyroresonance  heating 
within  the  ionosphere  would  involve  finding  the  sel 1-consistant  radiation 
pattern  of  the  heating  antenna  and  the  electron  distribution  function  by 
simultaneous  solution  ol  Maxwell  s  equations  and  the  Boltzmann  equation. 

The  complexities  of  this  problem  are  greatly  reduced  if  it  is  simplified 
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Figure  1.4,  Locations  and  orientation  of  antennas  with 
respect  to  the  airframe 


14 


to  one  dimension.  This  will  demonstrate  the  time-distance  dependence  of 
electron  temperature  upon  the  radiated  electric  field 

The  mode1  chosen  to  approximate  the  actual  situation  of  a  rocket-borne 
gyroresonance  heating  source  consists  of  a  plane  wave  incident  upon  an 
infinite  homogeneous  plasma.  The  direction  of  propagation  will  be  taken 
parallel  to  that  of  the  applied  magnetic  field,  and  it  will  be  assumed 
that  the  source  produces  circularly  polarized  waves  having  the  sense  of 
the  extraordinary  characteristic  waves. 

In  the  presence  of  a  gy rof requency  electric  field  of  extraordinary 
sense  of  polarization  the  aveage  energy  gained  by  an  electron  between 
collisions,  w,  is  equal  to  the  work  done  on  an  "average  electron” by  the 
field  in  the  same  interval  of  time  Based  upon  the  assumptions  of  the 
previous  section,  isotropic  collisional  scattering  implies  that  the  "average 
electron"  has  zero  momentum  immediately  following  a  collision  (i.e.  its 
energy  is ^thermal ized"  by  the  collision).  For  quasi-longitudinal  propaga¬ 
tion  of  the  gy rofrequency  field, 


2  2 

e  E  2 
w  =  <T*> 

2m 


(1.1) 


2 

where  E  is  the  field  intensity  and  ^  is  the  mean  square  free  time  for  an 

2  2 

electron.  Assuming  a  Poisson  distribution  of  free  time,  <t  >  =  2/v  and 


w 


2  2 
e  E 


mv 


(1.1a) 


(Derivations  of  ’.his  result  may  be  found  in  References  1  and  14.) 

The  average  electron  energy  in  the  presence  of  a  gy rof requency  field 
satisfies  the  equation 


t-  •  GMU  -  U  )  -  wv  (1  2) 

at  o 

where  U  =  average  electron  energy  (-  1/2  m  <v2>  =  3/2  kTe ) 

U  =  average  molecular  energy  (=  3/2  kT  ) 
o  o 

Letting  v  represent  rms  electron  velocity  and  i  r  i  v/v 

o  o’ 

(1.2)  (for  t  >  0  and  E  r  0  when  t  <  0)  is  given  by 


the  solution  of 
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V  C  _ 

O  /  >  2 

(v'  + 


v"2)  t  =  V'  ^  tanh  1  —■  -  tanh  -^y  )  -  v"  ^  tan  1 


,  V 

v  .  -1  o 

— rr  -  tan  —  rr 
V  V 


(1.3) 


where 


It  is  evident  that  v  increases  as  t  increases,  which  corresponds  to 
electron  temperature  with  time  during  the  gy rof requency  disturbance; 
we  call  this  a  heating  signal.  As  t - 


tan 


-1  v 


<1 


increasing 

thus 


for  all  v  and  therefore 


tanh 


-1  v_ 
v' 


->  oo 


or  v  must  approach  a  heated  equilibrium  value  v  .  Thus 

v  <  v  <a  v'  for  t  >  0 

o 

and  v - ^  v’  as  t- - ^  oo. 

We  now  make  the  assumption  that  electron  density  remains  constant  (i.e, 
neglect  effects  of  heating  on  diffusion,  ionization,  attachment  and  recombina¬ 
tion  rates).  The  heating  field  intensity  as  a  function  of  time  and  distance 
z  from  the  source  is 


E(z, t) 


E  exp 
o 


UI 

c 


,  t )  d£ 


(1.4) 
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where  E  is  the  field  at  the  source  and  "X  is  the  absorption  index  of  the  medium, 
o 

In  the  first  approximati on  %  may  be  assumed  to  be  that  found  from  the  steady 
state  Appel ton-Hartree  equations; 


%  = 


(1.5) 


where  «  is  the  plasma  frequency.  Now  from  Equation  Cl. 3)  and  Equation  (1.5) 
N 

we  know  %  as  a  function  of  time  only  when  E  is  constant,.  Therefore  it  is 
necessary  to  approximate  the  solution  to  Equation  (1.4)  in  the  following 
manner , 

It  is  assumed  that  the  source  of  the  disturbing  wave  is  turned  on  at 

t  =  0,  and  from  then  until  t  =  t  ,  an  interval  of  approximately  1/ Vq ,  it  is 

assumed  that  ^i(z,t)  =  %(z,0)  or  in  other  words  does  not  change  in  the 

interval  The  value  of  iE|  can  then  be  found  as  a  function  of  z  from 

Equations  (1  4)  and  (1.5)  and  the  increase  in  electron  velocity  at  various 

distances  from  the  source  (for  the  time  interval  from  t=Ot  t  =  t  )  can 

o  1 

be  found  from  Equation  (1.3).  From  these  data  and  Equation  (1.5)  a  curve 
of  X  versus  z  for  t  =  t  is  plotted  and  the  integration  of  Equation  (1.4) 
performed  graphically.  This  gives  |E(z,t  )l  which  is  held  constant  for  the 
next  interval  while  the  process  is  repeated  Successive  repetition  gives 
the  approximate  solution  to  Equation  (1  4)  for  later  values  of  time. 

Sample  calculations  of  this  type  have  been  made  for  altitudes  of  85  km 
and  92  km,  which  are  in  the  range  of  interest  for  the  rocket  experiments. 

The  assumed  parameters  of  the  ionosphere  for  these  altitudes  are  tabulated 
below 


Altitude 
85  km 
92  km 

The  values  of 
obtained  from 
Figures  1 
field  E  and  ol 


N 

7  0  x  10 

6 

8.4  x  10 


1  48  \  10 
4  20  \  10E 


T  °K 
o 

165.7 

168.4 


are  calculated  from  data  of  Waynick  ,  those  of  v  are 
N  o 

11  12 
Nicolet  ,  and  molecular  temperature  T  from  Minzner.  et  al 

o 

5  and  1  6  show  computed  curves  of  intensity  of  the  disturbing 


v/v  respectively  with  time  as  a  parameter  for  E  -  10  volts 
o  o 


-(VOLTS/METER) 
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Figure  i.o.  Field  strength  of  heading  wave  as  a  function  of  distance 

from  source  for  various  instants  of  time.  Field  strength 

at  source  E  =  10  vol ts/meter ;  altitude  =  85  km 
o 


0  100  Z-(METERS)  200  300 


Figure  1.6.  Relative  electron  velocity  as  a  function  of  distance 
from  source  of  heating  wave  for  various  instants  of 
time.  Field  strength  at  source  E  =  10  volts/meter; 
alt itude  =  85  km 
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per  meter  at  85  km,  Figure  1.7  and  1  8  show  E  and  v/v^  respectively  for  the 
same  time  intervals  and  at  the  same  altitude  but  for  E  =  1  volt  per  meter. 
In  Figure  1.9  and  1.10  E  and  v/v^  are  plotted  for  various  times  for  Eq  =  1 
volt  per  meter  at  an  altitude  of  92  km. 

1.6  Relaxation  of  Heated  Electrons 

At  the  end  of  a  gyroresonant  heating  pulse  it  is  presumed  that  the 
medium  has  reached  a  heated  equilibrium  state.  At  this  time  the  electron 
rms  velocity  is  roughly  v'  given  in  the  preceeding  section  (where  E  is  the 
electric  field  pattern)  Equation  (1.2)  now  becomes 


dtr 

dt 


Gv  ( U 


U  ) 
o 


0 


(1.6) 


which  has  the  solution  (when  the  substitution  U  =  1/2  mv  is  made) 


v 

v 

o 


Gv  t 

v  .  °  . 

—  ctnh  — —  +  1 

v  2 

o _ 

Gv  t 

v'  o 

—  ♦  ctnh  — — 

v  2 

o 


(1.7) 


This  indicates  that  1 /Gv  is  the  relaxation  time  constant  for  the  medium. 

o 

1.7  The  Gyroresonance  Heating  S ignaj_ 

The  criteria  for  choosing  pulse  width,  pulse  repetition  frequency, 
transmitter  power  and  antenna  system  for  the  gyroresonance  heating  signal 
are  as  follows.  It  is  desired  that  the  heating  spread  over  as  large  a 
volume  as  possible  In  addition  the  time  interval  between  pulses  should  be 
long  enough  to  permit  the  medium  to  relax,  this  insures  that  the  altitude 
variation  of  electron  density  and  collision  frequency  are  valid  experimental 
parameters . 

The  results  of  Section  1.5  indicate  that  heating  at  large  distances 
from  a  gyrof requency  source  occurs  much  more  slowly  than  near  the  source. 

The  heating  pulse  width  for  the  initial  experiment  was  chosen  to  be  500 
microseconds.  This  was  determined  mainly  from  the  practical  considerations 
involved  in  designing  equipment. 

As  pointed  out  in  Section  1.6  the  time  constant  for  the  decay  of 
electron  energy  in  1/Gv^  Therefore  the  pulse  repetition  frequency  should 


Figure  1,7.  Field  strength  of  heating  wave  as  a 

function  oi  distance  from  source  for 

various  instants  of  time.  Field  strength 

at  source  E  =1  volt/meter;  altitude  =  85  km 
o  ’ 
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Figure  L.H.  Relative  electron  velocity  as  a  function 
ol  distance  lnim  source  of  heating  wave 
lor  various  installs  of  time.  Field 
strength  at  source  K  =  1  volt/meter; 
altitude  -T  Kf)  km 
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Figure  1.10.  Relative  electron  velocity  as  a  function  of  distance 
from  source  ol  heat ing  wave  for  various  instants 
of  time.  Field  strength  at  source  F.  =  1  volt/meter; 
alt  it  tide  =  02  km  ° 
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be  less  than  Gv  At  100  km  v  10  and  for  G  -  10  ,  we  have  the  result 

o  o’ 

that  the  repetition  frequency  shiuld  be  less  than  100  pulses  per  second. 

From  system  synchroniza' ion  considerations  a  frequency  of  75  pps  was  chosen 
with  a  45  pps  free  running  capabili'y  should  the  synchronization  fail  This 
also  would  permit  analysis  ot  data  from  successive  pulses  as  functions  of 
geomagnetic  field  aspect  of  the  bea'ing  an'enna  based  on  a  rocket  rotation 
rate"  of  2  rps. 

The  heating  transmitter  output  power  level  is  an  arbitrary  parameter 
in  that  the  antenna  radiating  efficiency  depends  upon  its  impedances  which 
may  be  expected  to  be  a  function  of  al’itude  magnetic  aspect  and  degree 
of  electron  heating  Therefore  the  switching  of  the  output  power  through 
several  levels  seems  the  prudent  approach  For  the  initial  experiment  a 
transmitter  with  400  watts  peak  output  power  was  designed  A  step  attenuator 
between  the  transmitter  and  the  an'enna  provided  power  lpvels  of  4,  40  and 
400  watts.  Each  level  was  maintained  for  45  pulses  (or  ,6  sec)  which  would 
be  sufficient  time  for  one  roske'  rota' ion 

The  problems  involved  in  d*  signing  a  heating  antenna  system  do  not 
appear  to  have  corollaries  in  fre*  space  The  medium  is  anisotropic  and 
during  heating  inhomogeneous  It  is  furthermore  dependent  upon  parameters 
which  are  functions  of  alti'ude  The  antenna  system  design  for  the  initial 
experiment  consisted  of  a  pair  of  base  loaded  fiberglass  whip  antennas  each 
10  feet  long  The  impedance  of  shirt  antennas  in  an  anisotropic  medium  is 
discussed  in  the  work  of  K  Balmain  m  Chapter  2  of  this  report,  and  the 
heating  antenna  system  is  discussed  in  Quarterly  Progress  Report  No.  10. 
Briefly  the  effect  of  the  i  tno.spht  n  a'  90  km  is  ' o  cause  a  resistance  of 
the  order  of  2000  ohms  to  appear  in  shun'  with  the  antenna  terminals  after 
electron  temperature  has  been  increased  bv  a  factor  of  6  However  the  effect 
of  increased  electron  density  dm*  t  >  i  mizu'ion  would  be  to  lower  this 
resistance  Therefore  a  cumpr  omi value  of  1000  ohms  attributable  to  gyro- 
resonance  heating  losses  was  isn.ni  and  the  an'enna  matching  network 
designed  accordingly 

A  block  diagram  oi  'he  heating  'tansn-i'ier  subsystem  is  shown  in 
Figure  1  11.  The  synchrnnizi'  or.  >f  ’  h**  various  functions  is  discussed  in 
Section  110 
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Figure  1.11.  Heating  transmitter  subsystem  block  diagram 
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Data  from  which  antenna  impedance  and  radiated  power  could  be  computed 
was  obtained  by  measuring  four  RF  voltages  Three  (v^,  and  v^)  were 
along  a  50  ohm  coaxial  cable  at  one  eighth  wavelength  intervals  and  the 
fourth  (v^)  was  the  antenna  voltage.  These  probes  were  sampled  successively 
in  such  a  manner  that  for  each  heating  pulse  a  different  voltage  was  measured. 
A  gain  control  actuated  by  the  position  of  the  step  attenuator  provided  a 
constant  range  of  data  level. 

The  heating  signal  electric  field  data  acquisition  consisted  of  the 
measurement  of  the  RF  voltage  between  two  parallel  plates  near  the  upper 
tip  of  the  nose  cone.  This  configuration  was  chosen  because  its  impedance 
was  affected  less  than  a  dipole  by  variation  of  the  dielectric  properties 
of  the  medium.  Again  the  gain  of  the  amplifier  was  adjusted  in  accordance 
with  step  attenuator  position. 

1.8  Detection  of  Electron  Heating 

In  "Luxembourg  effect"  experiments  the  electron  heating  is  commonly 
detected  by  observing  the  modulation  produced  on  a  relatively  weak  signal 
traversing  the  disturbed  portion  of  the  ionosphere  The  modulation  is  an 
integrated  effect  of  the  absorption  over  the  entire  path  and  has  the  form 


Modul at l on 


T£tz,  i) )  dz 


(18) 


where  'Y-  (z)  is  the  undisturbed  absorption  index  and  V(z,t)  the  disturbed 
o  ’ 

absorption  index. 

For  a  sensing  wave  of  ordinary  polarization  and  propagating  parallel 
with  the  magnetic  field  in  a  homogeneous  plasma 


*<z, t ) 


1m 


X _ 

-  l  Z  ( z.  t ) 


(19) 


2 

where  X=  (u>  /^) 

Y  =  u.  /u,  -eB/m~) 

c 

Z(z  t)  =  I  v  /u,Tv(z,t)/v  1 
’  'o'1  ’  o 1 

The  dependence  of  on  v/v  is  responsible  for  the  cro^s  modulation  as  was 

V 

shown  by  Bailey  and  Murtyn  A  plo-  of  u.  •  c  TC  versus  v'v  or  T  a  i  85  km 

o  e 

for  a  sensing  frequency  1  5  times  the  gy rofrequencv  is  shown  in  Figure  1  12 


UJ/c  X  NEPERS  per  METER 
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v/v0 


w/c  1C  as  a  function  of  v/v  or  T  at  85  km  for  sensing 
wave  with  w  =  1 . 5  (=  2,1  me) 


Figure  1.12. 
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Note  that  the  absorption  coefficient  has  a  maximum  value  this  occurs  when 

Z/@L  +  Y)  =  1.  Therefore  the  maximum  absorption  at  lower  altitudes  will  occur 

for  smaller  value  of  v/v  and  for  higher  altitudes  at  greater  values  of 

o 

v/v  . 
o 

The  sensing  wave  frequency  and  time  of  day  for  the  experiment  were 
decided  on  the  basis  of  calculations  of  the  attenuation  and  derivative  of 
attenuation  with  collision  frequency  for  characteristic  waves  traveling 
parallel  to  the  magnetic  field  This  work  is  reported  in  Quarterly  Progress 
Report  Nos.  5  and  6  and  will  be  briefly  summarized  here 

The  night  experiment  was  ruled  out  because  the  ordinary  wave  is  not 
attenuated,  nor  is  it  sufficiently  affected  by  changes  of  collision  fre¬ 
quency  by  a  factor  of  10.  (as  would  be  expected  during  heating) 

The  daytime  calculations,  using  the  electron  density  and  collision  fre¬ 
quency  models  shown  in  Figure  1  lb,  indicated  that  the  sensing  wave  should 
be  near  2  me.  A  plot  of  the  expected  total  attenuation  for  the  ordinary 
wave  as  a  function  of  altitude  is  shown  m  Figure  1  13  The  extraordinary 
wave  is  attenuated  more  than  90  db  at  80  km 

A  broadcast  transmitter  was  given  to  i he  University  of  Illinois  by  the 
Collins  Radio  Company  for  use  in  these  experiments  as  a  sensing  wave  trans¬ 
mitter,  Having  been  modified,  it  is  now  capable  of  continuous  13  kw  output 
at  2.02  me  (the  assigned  sensing  frequency)  and  is  installed  m  a  semi-trai'ler 
At  the  rocket  it  is  desirable  that  the  sensing  antenna  be  unaf¬ 
fected  by  va-riAtions  of  the  dielectric  properties  of  the  medium.  The 
structure  which  gives  the  closest  approximation  to  this  requirement  is  a 
small  loop  antenna  This  is  also  aerodynamical ly  attractive  since  it  can 
be  mounted  within  a  fiberglass-epoxy  ogive  One  would  expect  the  efficiency 
of  such  an  antenna  to  be  low.  In  free  space  the  Q  of  a  lossless  circular 
loop  is  given  by 


Q 

i  > 


where  b  is  the  radius  of  the  1  •>  >p  and  a  the  taditts  of  t  Iil  conductor  For  2  me, 
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letting  b  =  2m  and  a  =  003  m  it  is  fuund  that 

Q  “  1 .4  x  107 
o 


However,  to  be  able  to  observe  details  in  the  sensing  wave  modulation 
during  the  heating  pulse  (pulse  duration  of  500  microseconds)  a  receiver 
output  bandwidth  of  DC  to  12  5  kc  was  chosen  (This  is  compatible  with  the 
IRIG  standard  70  kc  -  40 $  telemetry  subcarrier  channel  )  The  output  fre¬ 
quency  range  neccessi fates  an  RF  input  bandwidth  of  50  kc  or  the  antenna 
must  have  a  loaded  at  2  me  of 


The  resulting  antenna  efficiency  is  just  Q^/Q  ^3  x  10 

From  Figure  1.13  we  see  that  the  ordinary  sensing  wave  is  attenuated 
27  db  in  propagating  from  the  ground  to  100  km  There  is  also  a  loss 
of  3  db  in  the  extraordinary  wave  and  a  loss  of  3  db  due  to  the  linear 
polarization  of  the  receiving  loop  antenna.  If  the  sensing  transmitter 
output  power  were  10  kw  the  33  db  loss  would  make  it  equivalent  to  a  5  watt 
transmitter  in  free  space  with  identically  polarized  antennas  The  power 
density  expected  at  100  km  would  be  4  ■  10  11  watts  per  square  meter  the 
power  delivered  to  the  receiver  input  is  given  by 


P 


o 


GX 

4  a 


\ 


(Power  density) 


AV 


3  6  x  10 


wa  1 1  s 


For  a  receiver  input  impedance  of  100  ohms  the  resulting  voltage  at  the 
receiver  input  is  approximately  6  microvolts  This  is  adequate  for  detection 
with  a  high  signal  to  noise  ratio  using  stale  of  'he  art  techniques.  The 
sensing  wave  receiver  used  in  the  initial  gy i oresonance  experiment  was 
developed  for  this  purpose  bv  Kellogg  Division  of  I T  &  T 

The  rocket  rotation,  sensing  receiver  output  and  its  AGC  voltage  are 
sufficient  data  from  which  the  polarization  ellipse  of  the  sensing  wave  may 
be  determined  as  a  function  of  altitude  in  the  D  region  This  gives  the 
ratio  of  ordinary  to  extraordinary  wave  intensity  which  is  indicative  of 
•-the  properties  ol  the  1  >wt  t  region  of  i  h<  l  onosphet  e 
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1.9  Telemetry  Data  Conditioning 

The  telemetry  channels  available  for  this  experiment 
70  kc  -  40%  subcarrier  oscillator  using  234  me  and  70  kc, 
(all  -  15%)  subcarrier  oscillators  at  225  me  These  were 
following  manner 


consisted  of  a 
40  kc,  and  22  kc 
allocated  in  the 


70  kc  -  40% 
70  kc  -  15% 
40  kc  -  15% 
22  kc  -  15% 


sensing  receiver  output 

heating  antenna  impedance  and  radiated  power 
heating  field  intensity  and  excitation  light 
sampled  data 


The  data  conditioning  for  each  of  these  channels  will  be  discussed 

(1)  70  kc  -  40%  On  the  basis  of  computations  of  cross  modulation 

contained  m  Quarterly  Progress  Report  No  7,  the  expected  varia¬ 
tion  of  sensing  receiver  output  due  to  heating  is  of  the  order  of 
40%  at  85  km  and  10%  at  92  km.  The  sensing  receiver  output  has  a 
frequency  response  of  1  db  fiom  DC  to  12  kc.  The  DC  component 
is  held  at  -  10  v  (-  1/2  v)  by  the  AGC.  (Thus  a  1  volt  variation 
corresponds  to  10%  change  in  the  RF  input  voltage.)  However  the 
0-5  v  limits  of  the  subcarrier  necessitated  some  data  conditioning. 
This  was  accomplished  by  biasing  the  receiver  output  with  a  6.7 
volt  battery,  resulting  in  a  telemetered  range  of  sensing  wave 
amplitude  fluctuations  of  -  17%  and  -33% 

(2)  70  kc  -  15%  The  heating  antenna  impedance  and  radiated  power 

data  involved  a  time  sharing  scheme  for  use  of  the  telemetry 
channel  More  specifically  the  measured  quantities  were  the  voltages 
at  three  points  along  the  transmission  line  feeding  the  antenna 
(separated  by  eighth  wavelengths)  and  the  antenna  voltage  The 
change*  of  altitude  and  magnetic  aspect  of  the  rocket  between 
successive  heating  pulses  is  small  Therefore  during  each  heating 
pulse  a  different  voltage  (and  pulse  shape)  was  measured  Four 
successive  pulses  would  provide  sufficient  data  to  compute  the 
variation  of  antenna  impedance  and  radiated  power  during  the  heating 
provided  the  other  simultaneous  measurements  (sensing  wave,  field 
intensity)  should  indicate  that  'he  four  heating  pulses  produced 

the  same  el  led. s 


29 


(3)  40  kc  -  15%  This  channel  telemetered  the  field  intensity  and 
excitation  light  in  a  time  sharing  mode  Time  sharing  was  accom¬ 
plished  by  a  relay  which  was  actuated  to  provide  measurement  of 
field  intensity  for  the  first  30  heating  pulses  and  of  excitation 
light  for  the  last  15  heating  pulses  for  each  power  level  (The 
heating  power  level  was  changed  after  every  45th  pulse. ) 

(4)  22  kc_^  15$  Data  which  were  expected  to  be  very  slowly  varying 

were  sampled  by  a  30  position,  5  rps  commutator.  Position  and 
corresponding  data  are  listed  below 

1  +5  v  frame  sync  pulse 

2  *5  v  frame  sync  pulse 

3.  heating  transmitter  subcommutated  data 
4  subcommutator  position 
5.  photomultiplier  tube  high  voltage 
6  x  magnetometer  bias  voltage 
7.  y  magnetometer  bias  voltage 
8  z  magnetometer  bias  voltage 
9.  x  magnetometer  output 
10,  y  magnetometer  output 

11  z  magnetometer  output 

12  sensing  receiver  output  DC  component 

13  sensing  receiver  AGC  voltage  expanded 

14  sensing  receiver  AGC  voltage 

15  fuel  valve  data 

16.  *  4  05  v  calibration  battery 

17.  heating  signal  step  attenuator  position 

18  heating  antenna  deployment  data 

19  heating  transmitter  peak  RF  voltage 

20  common  ground  reference 

21.  heating  antenna  peak  RF  voltage 

22  sensing  receiver  output 

23  sensing  receiver  output 

24  x  magnetometer  output 

25  y  magnetometer  output 
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26,  z  magnetometer  output 

27.  sensing  receiver  output  DC  component 

28  sensing  receiver  AGC  voltage  expanded 

29,  sensing  receiver  AGC  voltage 

30.  difference  frequency  receiver  integrated  output 

The  data  on  positions  3  and  4  was  obtained  from  a  stepping  switch  within 
the  heating  transmitter  package  It  was  advanced  one  step  with  each  revolution 
of  the  data  sampling  commutator  The  measurements  made  concerned  the  monitoring 
of  the  operation  of  the  heating  transmitter  These  were 

1,  -t-13  v  battery  (power  supplies) 

2,  +6  5  v  battery  (filaments) 

3.  7 00  v  power  supply 

4.  <-300  v  power  supply 

5  oscillator  current 

6,  driver  amplifier  current 

Provision  was  made  to  advance  this  ^wi'ch  and  to  measure  these  voltages 
through  the  umbilical  prior  to  launch.  This  was  to  insure  proper  operation 
of  the  heating  transmitter 
1 , 10  Synchronization 

The  data  acquisition  processes  which  have  been  described  required  a 
source  of  synchronization  for  the  entire  experiment  The  data  sampling 
commutator  was  one  deck  of  a  two  pole  switch  The  o'her  was  used  as  a  clock, 
which  consisted  of  a  5  rps  60  position  switch 

Among  other  things,  ’he  clock  advanced  a  12  position  stepping  motor 
which  drove  the  heating  signal  step  a'tenuaror  Each  revolution  of  the 
clock  advanced  the  stepping  motor  once, every  third  clock  revolution 
advanced  it  twice  (Three  revolutions  of  the  clock  stepped  ’he  motor  4  times  ) 
Switches  driven  by  the  s*ep  attenuator  motor  provided  logic  circuits  from 
which  other  synchronization  was  derived 

Figure  1  14  shows  (he  time  relati  >nships  >[  the  synchronized  functions 
of  the  experiment  for  one  complete  cycle  of  the  step  attenuator 

The  data  calibration  markers  were  impressed  upon  the  three  continuous 
data  channels  at  rhe  times  shown  Each  marker  consisted  of  3  cycles  of  a 
400  cps  square  wave  The  maximum  )f  '  h  i ->  signal  was  -4  03  v  tl  t  s  and  ’he 


Figure  1.14.  Time  relationships  of  the  synchronized  functions  of  the  experiment 
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minimum  0  volts  These  were  provided  for  calibration  purpose,  and  as 
indicators  of  the  heating  step  attenuator  position  Note  'ha'  one  marker 
occurred  just  prior  to  each  change  of  attenuation  and  that  a  second  marker 
appeared  when  the  -20  db  position  was  reached 
1.11  Experimental  Results 

The  initial  gy roresonance  probe  was  launched  from  Eglin  Air  Force  Base 
at  1318  hours  CST  on  1  May  1962  The  primary  objective  of  the  experiment, 
the  detection  of  gy roresonance  heating  of  electrons  in  the  ionosphere,  was 
successful.  However  the  data  sampling  switch  failed  a  few  milliseconds 
after  burnout  when  the  acceleration  reached  -1  4  G  (from  photo’ heodol l te 
dftta) „  The  stoppage  was  sudden  indicating  a  mechanical  failure,  probably  due 
to  foreign  matter  drifting  upward  into  the  gears  as  the  acceleration  became 
less  than  -1 

This  failure  resulted  in  the  loss  of  sampled  data  and  of  the  synchroniza¬ 
tion  functions  As  a  result  the  following  conditions  prevailed  after  burnout 
and  for  the  remainder  of  the  flight 

(1)  The  heating  transmitter  operated  a'  45  pulses  per  second 

(2)  The  step  attenuator  remained  in  the  0  db  position 

(3)  70  kc  -  40%  subcarrier  channel  functioned  normally  This  channel 
carried  the  sensing  wave  receiver  output 

(4)  70  kc  -  1 5%  subcarrier  channel  carried  no  data  The  data  switch 
failure  occurred  during  a  4  5  volt  reference  marker,  which  then 
persisted  throughout  the  experiment 

(5)  40  kc  -  15 %  subcarrier  channel  measured  the  heating  field  intensity 

(6)  22  kc  -  15%  subcarrier  channel  measured  the  sensing  receiver 
AGC  voltage 

A  composite  of  the  sensing  wave  modulation  at  various  altitudes  is  shown 
an  Figure  1,15  Local  equilibrium  values  of  plasma  radian  frequency  u  and 

-  >  hj> 

collision  frequency,  'from  References  11  and  13  respectively!  are  given 
The  distance  traveled  by  the  rocket  tie'^een  heating  pulse,  ls  alsi  indicated. 
Vertical  dashed  lines  represent  the  onset  of  hearing  pulse,  which  had  0  5 
millisecond  durations  The  time  between  pulses  *as  22  milliseconds  and 
therefore  this  scale  of  the  data  permit,  study  of  the  sensing  wave  only  during 
relaxation  Each  division  of  the  ampli-ude  scale  is  1  v  4'  which  corresponds 
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Figure  1.15a, 


Sensing  wave  receiver  output  showing  cross¬ 
modulation  at  various  altitudes 
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Figure  1.15c, 


Sensing  wave  receiver  output  showing  cross¬ 
modulation  at  various  altitudes 
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to  10^  variation  of  the  envelope  of  the  inpu'  Pf  1  rotr.  Us  me  ir-  <  ,lue  The 
data  is  biased  so  that  the  lower  limit  of  the  telemetry  (  r  jin  l  corresponds 
to  -33^  while  the'  upper  limit  represents  '.7'$'  nev:  n:o"  or  the  receiver  input 
It  is  evident  i rom  Figuic  1  15  that  except  lot  >etv  lew  altitude =  the 
relaxation  time  ol  22  milliseconds  was  no’  !■  ng  enough  Furthermore  the 
character  of  the  disturbance  «■  in  gene  ref  is  more  c  e>mp  1 1  c  at  e  d  then  would  be 
expected  if  electron  energy  decay  we  t  e  the  elomi”ent  n  1  a\  i' ion  process 
involved  This  suggests  that  the  he  eting  powe  r  was  t e >< *  gie  it,  and  that  some 
of  the  cross  modulatior  mav  have  been  e  >uhi:  by  :o"./  ,*101  near  the  rocket 

There  are  apparent  difletenees  in  iscendmg  .no  te'eenni'g  cross  modula¬ 
tion  at  some  altitude's  The  heating  m'en-js  should  n  o  e  be  en  deployed  at 
T  e  76  seconds,  corresponding  to  62  I-  m  Ur  ,g  in  .v  n  1'  e  rre  1  ive  a  this  operation, 
resulting  in  1  different  interne  eei-l  igui  i’H'r  du'i-g  iHui'  ana  descent  , 

At  higher  altitudes  the  dille-tenees  in.  stri-.ng  *  ,\c  ,pe  s  mav  haic 
been  due  to  rocket  attitude  lh  ,t.  ;s  the  he  .tmg  i-m*-.',  v  id:  a*  ion  patterns, 
and  hence  the  spatial  electron  <n<-’gv  distribution  rune'io-s  wi’h  respect 
to  the  sensing  wave  propagation  vee'or  would  p  >  eb  if .  v  "e<’  be  ’re  same  lor 
ascent  and  descent  The  elite*  ol  mtgneiie  I'Ci  e  1  upon  1  gv»r  !  *t  quenty 
antenna  radiation  pti’etn  nuv  be  found  m  'he  work  t  It  V.tti  ,  .'ieludea 
in  Quartetly  Piogress  Fepe  rt  \<> 

Some  int  e  re  st  ing  t  e  -u  1’  s  we  r  e  eil'i  el  ne  d  item  *  he  D'  .  ill  /  it  in"  ellipse 

measurements  for  1  he  sensing  wave  ,  particularly  'uring  a-een*  wile  r  e  the 

rocket  attitude  was  liitlv  stable  Ire  it. eke'  te'i'iem  t,te  »'a<  roughlv 
one  revolution  in  8  suwias  p  •.  63  k  the  >e:  •:"t;  ■  'c  po .  r:r.a'ion 

was  approximately  linen  Hn'ive  1  ,  et  6"  Im  11  h  <n  bitotne  ilm-t  circuit! 

This  indicates  that  the  electron  eien-itv  ne  a*  65  km  1  so*  negligible  .  ind 
is  substantiated  by  the  f  ie  '  tha*  ero  -  modu  1  e'  ion  be  g  ,n  t()  ,pp,  >  62  km 

A  quant  1 1  at  1  ve  ini  1  V’s  1  -  n!  the  e  1  e  e  •  »  on  de  ns  1  t  v  a’ld  e  ■  ,  '  1  -  1  >>’'  !  teque  ne  V 

necessary  to  produce  'he  ob-cived  e  h  mge  s  m  'he  n  el  •  .’it  ol  ’*e  ellipse- 

will  be  made  al  1 1  1  iurthi  '  u  at  1  1  e  due  *  1  on  <>:  ’he  'tit  trn  •  t  » t  1  - 1  ”s  ;  -g  w  ,ve 

receiver  AOC  volt  agt  is  pit  tome  t'  by  APt-t 

During  isct  nt  i"  'ne  -1 8  1  1  t  o  6<l  k  O'  '  t  g  1  o«  «.  hi  "  1  •  1  1  '  hi  rut  it  .in  t.| 

the  rocket  c  ails  1  et  1  hi  limn  p>  >  1  /  •  1 1  •’  -  1  '  1  >  1  e  e  .  •  .  "g  c  •  1  • 11  a  .n  d 

the  sensing  wevi  to  pi~s  irjoagh  c  or'-  -g*-»  e  .  1  •  :  m- 
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the  antenna  polarization  pattern)  an  unusual  form  oi  modulation  occurred. 

Figure  1,16  shows  the  vicinity  of  the  null  at  57  km  and  is  typical  of  all 
of  the  nulls  in  this  range  of  altitudes  .  (The  general  shape  of  the  receiver 
output  is  due  to  the  AGO  which  had  a  time  constant  of  about  0  2  second) 

The  modulation  is  characterized  by  negative  going  pulses  coincident  with 
the  heating  pulses  which  become  apparent  25°  of  rotation  prior  to  the  null 
and  increase  in  amplitude  as  the  null  is  approached.  Within  one  or  two  de¬ 
grees  of  the  null  the  sense  of  the  pulses  change  to  positive.  The  modulation 

o  o 

becomes  negligable  25  past  the  null.  In  the  remaining  130  between  nulls 

there  is  no  modulation. 

This  suggests  an  explanation  based  upon  the  modulation  of  the  direction 
of  polarization  of  the  sensing  wave,  similar  to  Faraday  rotation.  For  example, 
if  each  heating  pulse  rotated  the  sensing  wave  momentarily  in  a  direction 
opposite  to  that  of  the  rocket^  the-  receiver  input  would  momentarily  decrease 
prior  to  the  null  and  increase  after.  Assuming  that  the  antenna  polariza¬ 
tion  pattern  is  sinusoidal  the  percent  deviation  of  the  receiver  input  would 
be 


Deviation  100“*  tan  0  sin  6q,  (1  10) 

where  0  is  the  angle  between  the  null  of  the  receiving  antc-nm  and  the 

polarization  of  the  sensing  wave,  and  6cf,  is  the  change  of  angle  of  the  sensing 

wave  polarization  due  to  a  heating  pulse.  Note  that  a  very  small  value  of 

modulated  rotation  ( )  can  produce  a  large  and  easily  detectable  deviation 

near  the  null  (0  0)  The  observed  diti  fits  equation  (1  10)  fairly  well; 

c  -  o 

and  the  case  shown  in  Figure  1  16  corresponds  to  uct>  :  1 

The  angle  of  rotation  of  the  sensing  wave  at  altitude  h  measured  clock¬ 
wise  as  viewed  from  the  ground  is  given  by 

h 

<b  7T  \  (H  -  H  )  dz  (1.11) 

P  2c  J  e  o 

o 

where  h-  and  are  respectively  the  real  parts  oi  +  he  qu  isi-longitudinal 
refractive  indicics  of  the  ixtraornirtry  ind  ordin  irv  eh  tr  id  eristic  waves  . 
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In  the  altitude  range  of  48  km  to  60  km  collision  frequency  is  large  while 
electron  density  is  usually  assumed  to  be  small  .  Under  these  conditions 
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oj  to 


(1 .12) 


Magnetic  aspect  data  during  the  initial  part  of  the  flight  show  that 
the  rocket  rotation  was  clockwise  when  viewed  from  the  ground.  Therefore 
the  polarization  rotations  caused  by  the  heating  pulses  were  counterclock¬ 
wise,  or  momentary  decreases  in  the  normal  Faraday  rotation  of  the  region. 
Prom  Equation  (1.12)  this  would  be  expected  if  v  were  increased^ due  to  in¬ 
creased  electron  energy  produced  by  the  heating  pulses 

The  ionospheric  conditions  in  the  48  km  to  60  km  region  necessary  to 
explain  the  observed  sensing  wave  rotation,  if  the  above  assumptions  are 
correct, are  somewhat  different  from  generally  assumed  properties,  For  in¬ 
stance,  it  is  necessary  that  the  electron  density  be  from  100  to  400  elec¬ 
trons  per  cubic  centimeter,  which  is  one  order  of  magnitude  high.  In  ad- 

-4 

dition,  the  fraction  of  energy  lost  per  collision  2  x  10  ,  which  is 

almost  an  order  of  magnitude  low  It  is  also  required  that  v  2  x  10? 

at  60  km  and  that  the  radiated  heating  power  be  approximately  100  watts. 

The  latter  condition  is  questionable  unltss  the  whip  antennas  were  deployed 

prematurely.  The  sudden  attenuation  of  the  extraordinary  component  of  the 

sensing  wave  at  65  km  would  have  re  quire  a  the  same  order  of  magnitude  of 

7 

electron  density  and  a  collision  frequency  of  less  than  10  Therefore 
observations  of  modulated  wave  rotation  and  of  extraordinary  wave  attenua¬ 
tion  both  lead  to  the  same  ionosphe-*ri c  properties  near  60  km 
1.12  Conclusions 

The  preliminary  analysis  of  the  telemetered  data  obtained  1  rc-m  the 
initial  gyroresonancc  probe  shows  that  the  prime  objectives  of  the  experi¬ 
ment  were  accomplished  Cross  modulation  of  a  sensing  wave  due-  to  electron 
heating  in  the-  E  layer  was  measured  In  addition  it  was  demonstrated  that 
there  is  an  anisotropy  in  the  45  km  to  70  km  region  ol  the  ionosphere  which 
can  be  detected  by  two  very  sensitive  me  t heeds  the  change  oi  axial  ratio 

of  the  sensing  wave  pohnza'iein  ellipse  with  altitude,  ind  the  modulation 
of  the-  rotation  of  the  sensing  wive  with  electron  heeling 
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It  should  be  noted  that  these  are  new  methods  for  probing  the  ionosphere 
and  may  be  extended  to  other  altitude  ranges  by  suitable  ehoice  of  sensing 
frequencies.  The  following  experiments  ,'i  proposed  to  measure  the  proper¬ 
ties  of  the  45  to  75  km  region 

(1)  Polarization  ellipse  me  asure  me  "  t  TMs  experiment  would  involve 
the  simultaneous  detection  ol  the  polarization  ellipses  of  four 
sensing  waves  in  the  1  8  to  3  me  range  This  would  give  over¬ 
lapping  coverage  ol  the  55  to  75  km  region  and  provider  both  elec¬ 
tron  density  and  collision  frequeney  data  rF’G  channels  f  and  E 
would  be  required  for  telemetry  with  two  receivers  multiplexed  on 
each  channel 

(2)  Modulation  of  sensing  wave  rotation  in  this  experiment  c-lc-ctron 
heating  would  be  accomplished  bv  use  of  a  pulsed  2  me,  400  kw 
transmitter  located  on  the  ground  (The  10  kw  sensing  wave  trans¬ 
mitter  used  with  the  gyror c son me c  experiment  can  be  modified  for 
this  purpose  )  By  heating  in  rhis  manre t  i  uniform  5*  modulation 
of  electron  temperature  ,  m<a  'hi  re  lore  of  wave  rotation  should  be 
obtained.  A  sensing  wave  would  be  used  te  cctect  this  rotation 
just  ac  in  the  gvrorc^onarue  experiment  ana  the  receiver  output 
would  be  telemetered  through  a  70  kc  _  40^  charnel 

The  study  of  the  results  of  the  fir-t  gvrorc so" ante  experiment  is  in¬ 
complete,  partially  because  data  reduction  bv  AP  (  h  net  been  finished 
In  addition  it  is  desire  able  that  these  results  be  correlated  with  those  of 
a  second  expc-iimont  Much  has  bee"  u  rmed  m  the  process  ol  data  anilysis 
which  will  lead  to  impi  ovc  me  nt  s  in  the  design  of  the  '-eeond  gvroresonance  probe 
For  example  it  is  evident  th  it  the  radiated  heating  power  should  be  reduced, 
and  that  the  time  between  pulses  ba  ina  u  ,std  to  permit  a  more  complete  re¬ 
laxation  Also  a  probe  ♦(>  measure  local  electron  density  ana  collision 
frequency  and  their  relaxation  rates  between  pulses  ^nould  be  included. 
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2  THE  SHORT  DIPOLE  ANTENNA  AS  AN  IONOSPHERIC  PROBE 

K.  G.  Balmain 


2^1  Introducti on 

2,1,1  Purpose  and  Description  of  Research 

The  purpose  of  the  research  covered  by  this  part  of  the  report  is  to  study 
the  characteristics  of  a  short  electric  dipole  antenna  mounted  on  an  ionospheric 
rocket  and  to  provide  information  helpful  in  interpreting  the  results  of 
experiments  using  such  an  antenna.  Since  only  the  final  phase  of  such  a  study 
can  be  carried  out  in  the  ionosphere,  it  was  decided  to  perform  a  series  of 
laboratory  experiments  to  measure  the  impedance  of  a  short  probe  in  an 
anisotropic  plasma, 

A  program  of  experimental  and  theoretical  research  has  been  carried  out 
with  the  result  that  the  near  fields  and  impedance  of  a  short  electric  dipole 
in  an  anisotropic  plasma  can  be  predicted  with  fair  accuracy  at  low  signal 
levels  and  estimated  at  high  levels.  The  theoretical  work  makes  use  of  the 
quasi-static  approach  in  which  it  is  assumed  that  the  AC  electric  field  near 
the  antenna  is  much  stronger  than  the  AC  magnetic  field  The  theory  has  been 
checked  by  laboratory  experiments  which  were  carried  out  in  the  afterglow 
plasma  of  a  pulsed  DC  discharge  in  helium  and  neon  gases,  In  the  experiments 
the  antenna  was  a  short  monopole  probe  inserted  in  the  anode  end  of  a  cylindri¬ 
cal  glass  discharge  tube  The  plasma  was  made  anisotropic  by  a  DC  magnetic 
field  parallel  to  the  monopole  axis  Measurements  of  the  probe  impedance 
agreed  substantially  with  calculations  using  the  quasi-static  theory 

The  quasi-static  theory  indicates  that  a  very  short  antenna  is  capable 

of  transmitting  appreciable  power  into  a  lossless  anisotropic  plasma.  In 

isotropic  media  the  quasi-static  theory  predicts  no  power  transmission  at  all 

because  a  magnetic  field  is  not  induced  by  an  irrotational  electric  field 

In  the  case  of  an  anisotropic  plasma,  however,  an  irrotational  electric  field 

does  induce  a  magnetic  field  and  thus  electromagnetic  phenomena  (such  as 

radiation)  can  occur  This  effect  is  the  same  as  that  studied  by  Trivelpiece 
1 5 

and  Gould  who  used  the  term  space  charge  waves  to  describe  it  (the  authors 
mentioned  used  only  the  source-free  equations)  The  condition  for  quasi¬ 
static  power  transmission  is  that  the  medium  must  be  hyperbolic  (in  the 
sense  that  the  differential  equation  used  is  hyperbolic  m  form)  Under 


t, 
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hyperbolic  conditions  the  input  impedance  of  the  short  antenna  has  a  resistive 

component  and  the  near  fields  are  concentrated  close  to  a  conical  surface  with 

apex  at  the  antenna  and  axis  parallel  to  the  magnetic  field.  It  should  be 

noted  that  quasi-static  theory  has  been  used  extensively  in  the  study  of 

ferrites  where  the  analogous  propagation  modes  are  referred  to  as  "magneto- 

16  17 

static  modes"  (for  instance,  see  Trivelpiece  et  al .  or  Joseph  and  Schlomann  ). 

The  major  experimental  problem  is  electron  density  non-uniformity  which 

influences  the  antenna  impedance  considerably.  Ambipolar  diffusion  to  the 

antenna  surface  makes  the  plasma  non-uniform  throughout  its  volume  and, 

immediately  adjacent  to  the  surface,  there  is  formed  an  ion  sheath  in  which 

the  electron  density  drops  sharply.  In  the  experiments,  application  of 

positive  DC  bias  was  helpful  in  minimizing  this  effect  by  collapsing  the 

sheath.  In  accordance  with  ion  sheath  theory  (Langmuir  probe  theory)  the 

sheath  thickness  was  increased  by  RF  heating.  With  the  sheath  collapsed  RF 

heating  produced  a  measurable  change  in  input  impedance  when  the  input  power 
-4 

was  as  low  as  10  watts  Despite  the  collapse  of  the  ion  sheath  with  bias 
and  the  use  of  low  RF  input  power  to  avoid  heating,  the  plasma  close  to  the 
antenna  was  believed  to  be  non-uniform  due  to  ambipolar  diffusion.  It  is  shown 
analytically  that  a  non-uniform  electron  density  tends  to  magnify  the  apparent 
collision  frequency  (as  measured  by  a  short  RF  probe) 

It  should  be  emphasized  that  many  of  the  experimental  problems  encountered 
were  not  anticipated  and  consequently  the  apparatus  had  to  be  modified  several 
times.  The  vacuum  system  was  rebuilt  to  make  pumping  and  gas-changing  easier 
and  to  reduce  contamination.  The  electrode  design  was  modified  to  permit  DC 
biasing  of  the  antenna  probe  with  respect  to  the  plasma  and  to  make  the  pulsed 
discharge  more  dependable.  In  addition,  the  discharge  tube  was  made  smaller 
in  terms  of  wavelengths  when  it  was  decided  to  limit  the  research  to  short 
antennas . 

2. 12  Historical  Survey 

Rocket -mounted  probes  have  been  widely  used  to  measure  electron  density 
in  the  ionosphere  These  probe  experiments  can  be  classified  according  to 
the  type  of  applied  voltage  suggested  classifications  are  DC  experiments 
(Langmuir  probe),  low-level  RF  experiments  (which  do  not  perturb  the  plasma) 
and  high-level  RF  experiments  (which  do  perturb  the  plasma  locally)  A 
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number  of  DC  experiments  have  been  carried  out  by  a  group  at  the  University  of 
1 9-22 

Michigan  whose  work  centers  around  a  dumbbell-shaped  electrostatic  probe, 

23-25 

Low-level  RF  impedance  measurements  have  been  analyzed  by  Jackson  and  Kane 

26 

and  by  Pfister,  Ulwick  and  Vancour  who  considered  the  case  of  a  short  dipole 

in  an  isotropic  plasma.  In  addition  the  many  impedance  experiments  by  the 

University  of  Utah  deserve  special  mention  (see,  for  instance.  Haycock  and 
18 

Baker  )„  A  high-level  RF  approach  called  the  "Resonance  Probe"  method  has 

27  28 

been  described  by  Takayama  and  others  ’  who  make  use  of  an  anomaly  in  the 
DC  current  to  an  RF  probe.  This  anomaly  occurs  when  the  frequency  of  the 
applied  signal  passes  through  plasma  resonance  and  thus  it  permits  the  electron 
density  to  be  determined. 

The  theory  of  antenna  impedance  in  plasma  media  is  not  highly  developed. 

2 9  30 

Some  work  in  isotropic  media  has  been  done  by  King  and  others  ’  at  Harvard. 

31 

Katzin  and  Katzin  have  derived  impedance  expressions  for  a  dipole  of  arbitrary 
length  in  an  anisotropic  plasma  but  their  assumption  of  a  sinusoidal  current 


distribution  would  be  applicable  only  to  short  antennas.  A  quasi-static 

32 

approach  has  been  described  by  Kononov  et  al  who  considered  only  an  infinitesi¬ 
mal  dipole.  However  their  theory  appears  to  contain  an  error  which  prevents 
them  from  observing  the  resistive  part  of  the  antenna  impedance  mentioned 
earlier.  Other  investigators  ’  ’  have  used  quasi-static  theory  in 

their  studies  of  propagation  m  small,  anisotropic  structures  (both  plasma 
and  ferrite)  but  have  considered  only  the  source-free  solutions.  In  general 
one  may  say  that  the  impedance  theory  of  antennas  in  anisotropic  plasma  is 
in  need  of  further  development  and  that  laboratory  experiments  to  measure 
impedance  are  almost  non-existant , 

2.1.3  Summary  of  Report 

In  this  report  the  physical  factors  affecting  the  impedance  of  an  antenna 
are  discussed  with  special  attention  to  the  problems  peculiar  to  rocket  and 
laboratory  experiments.  The  quasi-static  theory  of  antenna  impedance  for 
anisotropic  media  is  developed  and  its  implications  discussed  The  laboratory 
experimental  apparatus  is  described  along  with  a  number  of  experimental 
results.  The  effect  of  a  non-uniform  electron  density  is  estimated.  The 
quasi-static  theory  is  applied  to  the  gyro- interact l on  problem  where  it  is 
used  to  estimate  the  impedance  of  the  rocket  antenna  and  its  field  intensity 
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as  measured  by  a  monitor  in  the  nose  cone.  Finally,  several  suggestions  for 
continued  research  are  put  forward 

An  appendix  contains  a  general  impedance  theorem  derived  by  Professor 
G,  A.  Deschamps  This  theorem  relates  the  impedance  of  an  antenna  in  a 
uniform,  isotropic,  conducting  medium  to  the  impedance  of  the  same  antenna 
in  free  space 

2.2  Physical  Factors  affecting  Impedance 
2. 2.1  Introduction 


The  fundamental  difficulty  in  predicting  the  impedance  of  an  antenna  in  a 
plasma  is  that  the  antenna  disturbs  the  surrounding  medium  Ions  and  electrons 
diffuse  to  the  surface  of  the  antenna  rendering  the  medium  inhomogeneous  and 
producing  a  positive  space-charge  sheath  (the  ion  sheath)  adjacent  to  the 
antenna.  An  RF  signal  applied  to  the  antenna  excites  the  medium  non-linearly 
and  heats  it  non-uni formly .  In  the  ionosphere  the  magnetic  field  of  the 
earth,  photo-emission  and  rocket  velocity  further  modify  the  antenna  environ¬ 
ment.  Since  these  factors  could  be  important,  they  will  be  reviewed  in  some 
detail.  This  chapter  is  largely  a  survey  of  published  work  together  with  some 
observations  made  in  this  laboratory 

2  2.2 _ The_I^ori  Sheath 

When  a  surface  is  introduced  into  a  plasma,  the  electrons  and  ions 
(thermally  agitated)  diffuse  to  the  surface  where  they  recombine  Since  the 
electrons  diffuse  more  rapidly  than  the  ions,  the  surface  acquires  a  negative 
charge  high  enough  to  establish  a  state  of  equilibrium  between  the  electron 
and  ion  currents  These  currents  are  equal  if  the  surface  i s  an  insulator 
or  a  "floating”  conductor  Thus  there  is  formed  a  region  of  positive  space 
charge  between  the  surface  and  the  adjacent  plasma  (which  is  electrically 
neutral).  We  shall  assume  the  boundary  between  the  plasma  and  the  ion  sheath 
to  be  fairly  well  defined  and  call  it  the  sheath  edge  ' 


The  ion  sheath  theory  of  a  plane  conducting  surface  will  be  presented 
this  is  usually  referred  to  as  the  plane  Langmuu  probe  (There  are  many 
references  to  Langmuir  probes  in  the  literature38,39,40  41 ’ 1 9< 20  '21 ’ 22, 34‘37  ) 
The  nomenclature  is  as  follows 


N  =  electron  and  ion  density  in  plasma 
-e  =  charge  of  electron 
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m  =  mass  of  electron 
M  =  mass  of  ion 
x  =  sheath  thickness 

V  =  voltage  across  sheath;  normalized  voltage  =  v  -  eV/kT^ 

€  =  1/3677  x  10^  =  8.85  X  10  ^2  farads/meter 

o 

k  =  Boltzmann's  constant 
T\  =  ion  temperature 
Tg  =  electron  temperature 
D  =  Debye-Huckel  shielding  distance 
JL  =  random  ion  current  density  crossing  sheath  edge 
J  =  random  electron  current  density  crossing  sheath  edge 
A  =  area  of  plane  probe 

For  a  Maxwellian  velocity  distribution,  the  random  current  densities  are 


(2.1) 


J. 

l 


(2.2) 


From  the  theory  of  the  space-charge-limited  plane  diode,  we  have 


J 

l 


(2  3) 


The  ions  which  cross  the  sheath  edge  all  reach  the  probe,  the  electrons,  how¬ 
ever,  are  repelled  by  the  negative  probe  potential.  Thus  the  currents  to 
the  probe  surface  are 


I 

e 


/kT  \1/2  -eV/kT 


3  \  277m 

)  • 
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fkM1/2 

A  Ne  I 
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\  m)  "2 

(2.4) 


(2  5) 


I . 

l 


A 


(2  6) 
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Let  us  assume  that  the  plasma  is  isothermal  (T  =  T  =  T)  and  introduce  the 

2ei/2  1 

Debye-Huckel  shielding  distance  D  =  (fkT/2Ne  )  which  is  the  distance  beyond 
which  a  plasma  is  effectively  shielded  from  the  field  of  a  point  charge 
Defining  v  -  eV/kT  and  using  Equations  (2  1)  and  <2  2),  we  may  write 
Equations  (24),  (2,5),  and  (26)  in  a  simpler  form 

1  =  A  J  e"V  (27) 

e  e 

I.  =  A  J  (2  8) 

i  l 


16JT1 ' 2  D2  3/2 

A  J.  T  1  " 

X 


(2.9) 


Equating  (2.8)  and  (2  9)  we  have  the  sheath  thickness  for  a  given  sheath 
voltage: 


2 

x 


1 67T1  '2  2  3/2 

-9-  D  V 


(2  10) 


If  there  is  no  net  current  from  the  probe  (or  if  we  are  considering  an 
insulating  surface)  we  may  set  I  -  1^,  giving 


(2.11  ) 


Substituting  Equation  (2  11)  in  Equation  (2  10)  gives  the  sheath  thickness 
for  an  insulated  probe 


x 


(2  12) 


The  shielding  distance  is  given  in  the  accompanying  graph  (Figure  2  1)  for 

T  =  290°  K  (for  which  D  -  D  ) 

o 

The  effect  of  bias  may  be  estimated  by  considering  Equation  i2  10) 
which  can  be  wri'ten 


1  '4 


x 


4tr 

3 


D 


(2  13) 
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where  V  is  the  probe-to-plasma  potential  If  the  probe  is  biased  negatively 

with  respect  to  a  large  electrode  immersed  elsewhere  in  the  plasma  (see 

Figure  2.2),  V  is  approximately  equal  to  the  bias  voltage  (the  difference 

being  the  sheath  potential  of  the  large  electrode,  usually  .1  to  .2  volts). 

Thus,  for  negative  bias  voltages  greater  than  about  one  volt,  the  sheath 

thickness  is  proportional  to  the  three-fourths  power  of  the  bias  voltage. 

For  an  insulated  probe,  Equation  (2  12)  states  that  x  is  proportional 
1/2 

to  T  as  long  as  T  =  T  =■  T  .  The  situation  is  not  so  simple  when  an 

l  e 

applied  RF  signal  selectively  heats  the  electrons,  leaving  the  ion  tempera¬ 
ture  unchanged.  For  this  case  (T  /  T^),  the  expression  for  sheath  thickness 
can  be  shown  to  be 


(2.14) 


which  indicates  that  x  is  approximately  proportional  to  1'  However, 

42  1/2^ 

Schulz  and  Brown  argue  that  J  is  proportional  to  rather  than  to 

1/2  1  e  1/2 
T\  .  If  this  were  the  case,  x  would  be  proportional  to  In  any 

case,  heating  (selective  or  not)  Lends  to  enlarge  the  sheath  This  was 

borne  out  by  early  laboratory  experiments  here  which  showed  that  RF  heating 

and  increased  negative  bias  produced  the  same  type  of  change  in  probe  impedance. 

One  may  obtain  an  estimate  of  sheath  thicknesses  to  be  found  in  practice 

o 

by  considering  the  case  T  -  290  K  for  which  v  =  40  V  The  shielding  dis¬ 
tance  for  this  case  is  D  and  is  given  in  the  accompanying  graph  Equation  (2.12) 

o 

becomes 


x  =  1  056  D  (in  --)  (2  15) 

For  the  case  of  ionized  helium,  this  becomes 

x  6  75  B  (2  16) 

o 

and  the  sheath  voltage  V  -  149  volts  Similarly  I  he  effect  of  bias  may  be 

estimated  by  re-writing  Equation  (2  Id)  as 

28  2  I)  V3A1 
(  > 


X 


(2  16a) 


51 


Since  the  theory  of  antennas  in  plasma  assumes  a  uniform  medium,  the 
use  of  positive  bias  to  collapse  the  sheath  should  be  considered  The 
probe  or  antenna  must  be  biased  with  respect  to  a  reference  electrode  with 
area  large  enough  to  permit  sheath  collapse  The  required  area  ratio  may 
be  estimated  as  follows  (refer  to  Figure  2  2) 

1  =  I  -  I  =  I  -  I  i2  17) 

ei  X1  *2  e2 

-V  -  l' 

A  ( J  e  1  -  J  )  =  A  ( J  -  J  e  2)  (2  18) 

1  e  l  2  l  e 


e 


(2,19) 


If  the  sheath  over  A  is  to  be  collapsed,  v  =  0  Also  J  /J  r  (m/M)  which 

1  1  l  e 

in  general  is  small  compared  to  unity  Thus  Equation  (2  19)  may  be  approxi¬ 
mated  by 


'2  20) 


For  this  equation  to  have  a  solution  it  is  required  that 


(2  21  ) 


For  example  consider  Helium  for  which 


Evidently  the  area  of  the  reference  electrode  should  be  much  greater  than 
the  area  of  the  biased  probe.  If  a  rocket  antenna  is  to  be  biased  with 
respect  to  the  rocket  body,  it  is  important  that  no  insulating  cover  be 
placed  around  the  anlenna  and  that  a  sufficient  area  of  the  rocket  body 
be  left  unpainted 
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40 

The  potential  distribution  in  a  plane  sheath  has  also  been  calculated 
Let  v  be  the  normalized  probe  potential  and  r)  the  normalized  potential  at 
any  point  z  in  the  sheath  (with  z  =  0  taken  to  be  the  probe  surface).  Then 
the  potential  distribution  is 

"0  =  4  tanh 

The  exponential  term  dominates  in  the  above  expression  and  emphasizes  the 
importance  of  the  shielding  distance  D 

So  far  only  the  plane  ion  sheath  has  been  considered  but  similar  analyses 
exist  for  cylindrical  and  spherical  sheaths.  The  basic  equations  are  sum¬ 
marized  in  Table  I.  Simplified  ion  current  formulas  are  given  for  the  "sheath- 
area-limited"  case  (when  all  the  ions  entering  the  sheath  are  collected  by 
the  probe)  and  for  the  'orbital-motion- 1 imited''  case  (when  the  space-charge 

in  the  sheath  does  not  influence  the  probe  current)  In  general,  graphical 

1  9  22 

methods  must  be  used  to  solve  the  equations  ’ 

2.2,3  Factors  affecting  the  Rocket  Antenna 

The  ion  sheath  already  discussed  is  modified  by  rocket  velocity,  photo¬ 
emission,  and  the  earth's  magnetic  field  These  effects  all  tend  to  make 
the  sheath  around  a  cylindrical  antenna  asymmetrical.  The  magnitudes  of 

the  effects  are  difficult  to  compute  precisely  but  can  be  estimated  An 

45 

excellent  review  on  this  subject  has  been  prepared  by  Zachary 

Let  us  compare  the  gyro- interact  1  on  rocket  vel oci ty v  with  the  rms 

electron  and  ion  velocities  (v  and  v  )  in  the  E  layer 

e  l 

5 

v  ^  1  2x10  m/sec 
e 

3 

v^06xl0  m'  sec  at  85  km 

3 

v  ^  0  6  x  10  m  se^ 

Since  vajvj  «  v^  it  appears  that  the  rocket  velocity  would  modify  the  sheath 

structure  somewhat  For  most  practical  purposes,  however,  it  is  believed  that 

rocket  velocity  effects  are  negligible 

The  magnitude  of  the  photo-emission  effect  is  by  no  means  settled 
22 

Boggess  estimates  that  the  photo-emission  current  density  in  the  E-layer 


-z/D  ,  v 
e  tanh  — 

4 


(2, .22) 


I 

i 

I 
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is  of  the  same  order  of  magnitude  as  the  random  thermal  electron  current  density 

and  thus  photo-emission  could  have  a  considerable  effect  on  the  sheath.  How- 
45 

ever,  Zachary  reports  later  research  leading  to  the  conclusion  that  the  photo¬ 
emission  effect  is  negligible  up  to  an  altitude  of  1000  km, 

1 9 

Hok  et  al.  express  the  opinion  that  the  earth's  magnetic  field  might 
affect  Langmuir  probe  measurements  of  electron  density  but  not  measurements 
of  ion  density.  This  means  that  there  will  be  a  magnetic  field  effect  as 
long  as  there  is  appreciable  electron  collection  by  the  probe  (as  in  the  case 
of  an  insulated  probe).  In  the  gyro-interaction  experiment  the  antenna  will 
always  be  nearly  perpendicular  to  the  magnetic  field  and  thus  any  sheath 
effect  should  be  nearly  constant  as  the  rocket  spins. 

In  conclusion,  the  rocket  effects  mentioned  above  are  not  expected  to 
have  a  great  influence  on  the  antenna  ion  sheath  and  hence  on  its  impedance. 
Nevertheless  in  interpreting  the  rocket  flight  records  it  will  be  necessary 
to  watch  for  anomalous  results  which  can  be  correlated  with  velocity,  solar 
aspect  and  magnetic  field.  When  the  gyro-interaction  heating  antenna  is 
operating  at  a  high  RF  level  the  plasma  will  be  strongly  disturbed  and  the 
effects  discussed  should  be  reduced. 

2.2.4  Factors  affecting  Impedance  Measurement  in  the  Laboratory 
In  the  afterglow  of  a  pulsed  DC  discharge  the  time  sequence  of  events 
is  very  important  and  may  be  outlined  as  follows: 

1.  The  DC  pulse  terminates 

2.  The  ions  quickly  reach  thermal  equilibrium  with  the  gas  molecules 
(near  room  temperature) 

3.  The_  electrons  more  slowly  reach  thermal  equilibrium  with  the  gas 
molecules,  mainly  through  a  process  of  diffusion  to  the  container  walls. 

4.  If  the  electron  and  ion  densities  are  high  enough  the  dominant  plasma 
decay  process  becomes  volume  recombination.  Since  this  process  does 
not  involve  the  container  walls,  the  spatial  distribution  of  charged 
particles  tends  to  become  more  uniform  (the  diffusion  process  of 
step  3  has  left  the  distribution  highly  non-uniform) 

5.  With  the  charged  particle  density  much  decreased,  diffusion  again 
predominates  and  the  plasma  density  distribution  becomes  more  non- 
uniform. 

The  time  required  for  the  attainment  of  electron  thermal  equilibrium 
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(periods  3)  is  of  major  importance  and  has  been  studied  by  Dougal  and  Goldstein 
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For  helium  and  neon  this  time  constant  (1  )  is  given  as  follows  (for  gas 

e 

pressure  p  between  ,5  and  5.0  mm): 

Helium  t^  <  8.4/p  +  27  p,s  for  p  =  1  mm,  t  <  36  (is  (2  23) 

Neon.  t  <  150/p  +  90  us  for  p  =  5  mm,  t  <  120  us  (2,24) 

-  e  —  1  e  — 

Thus  in  time,  t  ,  after  the  discharge  pulse  termination  the  electrons  should 
be  near  room  temperature.  However,  the  persistence  of  metastable  excited 
states  may  extend  this  time  somewhat. 

When  recombination  dominates  (period  4)  the  plasma  is  most  nearly  uniform 
and  good  experimental  results  are  most  likely  to  be  obtained.  Nevertheless, 
some  ambipolar  diffusion  will  still  be  present  so  that  a  completely  uniform 
plasma  is  impossible  to  obtain.  Helium  and  neon  are  suitable  for  this  type 
of  experiment  because  their  collision  cross  sections  are  well  known  and 
nearly  constant.  However,  neon  is  to  be  preferred  because  of  its  higher 
recombination  coefficient,  lower  diffusion  coefficient,  and  lower  collision 
probability  (which  permits  operation  at  higher  pressures,  thus  reducing  con¬ 
tamination  problems).  To  summarize,  at  300°  K  and  1  mm  pressure  we  have 
(in  cgs  units) 


Coll  is  ion 

Recombi nat l on 

Ambipolar  Diffusion 

Probability  P 

Coefficient  a 

Coefficient  D 

Hel ium 

19 

1  7  \  10"8 

540 

Neon 

3.3 

2  1  \  10'7 

115 

where  a  ancj  d  are  as  given  by  Goldstein47  and  P  as  given  bv  Brown48  A  DC 
a  c 

magnetic  field  reduces  transverse  diffusion  but  does  not  affect  longitudinal 
diffusion.  Thus  when  the  magnetic  field  is  parallel  to  a  cylindrical  dipole, 
diffusion  to  its  surface  would  be  reduced  and  the  plasma  consequently  would 
become  more  uniform 

During  the  course  of  the  experimental  work  it  was  found  that  the  RF 
probe  surface  became  contaminated,  probably  because  the  probe  is  part  of 
the  DC  discharge  anode  A  yellowish  waxy  deposit  formed  on  the  probe  surface, 
giving  it  an  insulating  coat  which  disturbed  the  impedance  measurement  and 
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made  sheath  collapse  impossible.  The  deposit  probably  consists  of  oxides, 
sputtering  products  from  the  cathode,  and  hydrocarbons  from  a  small  amount 
of  vacuum  grease  used  to  seal  the  RF  probe  mount.  The  only  known  cure  for 
this  condition  is  frequent  cleaning  of  the  RF  probe. 

It  is  essential  in  impedance  measurements  to  understand  the  loss  pro¬ 
cesses  and  include  as  many  as  possible  in  the  theory.  It  is  easy  to  include 
any  loss  process  which  may  be  represented  by  an  equivalent  collision  fre¬ 
quency  v.  Electron-molecule  and  electron-ion  collisions  are  most  important 

in  this  work  and  we  may  write  v  =  v  +  v  .  The  appropriate  formula  areas 

em  el  46 

follows  (given,  for  instance,  by  Dougal  and  Goldstein  and  discussed  by 
49 

Pfister  ): 

v  =  —  v  P  p  (MKS  units) 

em  3  c  o 

3.62  X  106  N,  I  3.30  X  10®  T  3/2 

"ei  =  T  3/2  ^  I  N1/2 

e  L_  i  __J 


where  v  -*/ 8kT  /7Tm  =  average  velocity  of  electrons  with  Maxwellian  distribution 
e 

p  =  273/T  p  =  pressure  reduced  to  0  C 
o 

N.  =  ion  density 

l 

These  results  are  presented  graphically  in  Figure  2.3  (a  and  b) .  The  graph 

shows  that  electron-ion  collisions  must  be  taken  into  account  in  the  work 

12 

covered  by  this  report  since  the  electron  density  may  go  as  high  as  10  per 
cc.  For  a  convenient  reference,  a  graph  of  plasma  frequency  as  a  function  of 

electron  density  has  been  included  (Figure  2.4),  in  MKS  units,  the  relation 

2  V 

is  f*  =  80.5  N. 

N 

There  are  other  loss  processes,  however,  which  are  non-coll isional  in 

50 

nature  and  thus  their  effects  can  only  be  estimated.  Tonks  has  observed 
that  non-uniformity  broadens  any  plasma  resonance  effect  thus  indicating  apparent 
high  losses.  In  addition  the  non-uniform  plasma  in  a  metal-covered  glass  tube 
can  support  surface  wave  modes  which  carry  away  energy,  again  indicating 
apparent  high  loss.  The  problem  of  non-uniformity  of  the  electron  density  is 
considered  further  in  Section  2.5  in  an  attempt  to  interpret  the  experimental 
results. 
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2.3  The  Quasi-Static  Theory  for  Near  Fields  and  Impedance 

2.3.1  The  Necessary  Assumptions 

In  the  theory  to  follow  it  is  assumed  that  the  plasma  medium  around 

the  antenna  is  linear  and  uniform,  two  requirements  often  difficult  to  satisfy 

in  practice..  The  linearity  requirement  is  easy  to  satisfy  provided  that  the 

RF  input  voltage  can  be  made  less  than  the  average  thermal  energy  of  an 

o  27 

electron  (.04  volt  at  300  K) ;  this  is  out  of  the  question  for  a  practical 
transmitting  antenna  but  can  be  satisfied  for  the  case  of  a  receiving  antenna 
or  an  RF-energized  plasma  probe.  A  high  RF  level  not  only  renders  the  medium 
non-linear  but  also  raises  the  electron  temperature  non-uniformly.  In  a 
theoretical  analysis  these  effects  can  be  approximated  very  roughly  by  increas¬ 
ing  the  assumed  losses  in  th^  medium.  A  non-uniform  electron  density  arises 
from  diffusion  to  the  antenna  surface  (or  any  other  nearby  surface)  and  is 
very  difficult  to  deal  with.  The  ion  sheath  resulting  from  diffusion  close 
to  the  surface  often  can  be  collapsed  with  positive  bias,  however  this  tends 
to  raise  the  electron  temperature  Ambipolar  diffusion  cannot  be  controlled 
in  the  ionosphere  but  in  the  laboratory  the  type  of  gas  and  the  pressure  can 
be  selected  to  minimize  it.  To  conclude,  linearity  and  heating  are  not 
problems  at  low  RF  levels,  but  a  non-uniform  electron  density  cannot  be 
eliminated,  only  minimized 

2.3.2  The  Quasi-Static  Theory 

The  problem  of  calculating  the  near  field  of  a  dipole  antenna  (and  thus 
its  impedance.)  is  not  easy  when  the  antenna  length  and  the  properties  of  the 
medium  are  arbitrary.  The  problem  can  be  simplified  greatly  if  the  antenna 
is  much  shorter  than  a  wavelength  this  permits  the  assumption  that  the 
cumnt  taries  linearly  along  each  hall  of  the  antenng  vanishing  at  the  ends.  Furthermore, 
in  an  isotropic  medium  the  magnetic  field  is  negligible  compared  to  the 
electric  field  in  the  near  region  of  a  short  dipole.  This  means  that  the 
electric  field  can  be  expressed  as  the  gradient  of  a  scalar  potential  and 
thus  the  antenna  behaves  as  if  it  were  a  capacitor  If  this  quasi-static 
theoretical  approach  were  applicable  to  anisotropic  plasmas,  the  input 
impedance  of  a  short  dipole  could  be  computed  readily.  However,  the  fact 
that  the  DC  magnetic  field  forces  moving  electrons  into  circular  orbits 
suggests  that  the  dominant  electric  field  may  induce  a  magnetic  field  which 
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is  not  negligible,  thus  invalidating  the  quasi-static  theory. 

This  problem  may  be  clarified  by  considering  Maxwell's  equations 


V  X  E  =  -juo  ijl  H 

o 

V  X  H  =  jco  £  Re 

o 


where  the  permittivity  tensor  is 
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=  collision  frequency 


If  it  is  assumed  that  -jw  p^H  is  small  enough  to  be  negligible,  then  V  \  E  =  0. 
Consequently  E  =  -V  where  is  a  scalar  potential.  In  order  to  see  whether 
or  not  an  E  field  of  this  form  (irrotational)  induces  a  magnetic  field,  take 
the  curl  of  the  second  Maxwell  equation 

V  X  V  x  H  =  joo  €  V  X  R  E 
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W-H  -  =  -Jtu  €  V  X  K  v  t 

o 


Since  V  *  H  =  0,  this  equation  becomes 

V2^  =  ju^VxRVijj 


If  G  is  defined  as  G  =  jw  €  V  X  K  V  ih  this  becomes 

o 

V2H  =  G 


or,  if  the  vector  components  (rectangular)  are  considered  separately,  VH.  =  G^, 
"i"  indicating  x,  y,  or  z.  Now  IT  may  be  expressed  in  terms  of  G^  using  Greens 
theorem  (A  is  arbitrary). 


/, 


V2 A  -  A  V^H  .  ) 


dv' 


J(H.  V  A  -  A  V  H  )  •  n  da’ 

1 


If  the  surface  s  approaches  infinity  and  if  the  medium  is  slightly  lossy, 
the  right-hand  side  will  approach  zero.  If  A  is  specified  as  follows: 


;  T  -  r’  !  R 


then 


V2A  =  -47 T  6  ( r 


r’  ) 


Substitution  in  Green's  threorem  gives  the  components  of  the  induced  magnetic 
field 


H  (r) 

l 


1_ 

477 


G  (r'  ) 

l 

r  -  T'  ; 


v 


dv 


Since  G  is  non-zero  in  general,  there  exists  an  induced  magnetic  field  H  which 
must  be  small  if  the  quasi -si at ic  theory  is  to  be  valid, 
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If  the  medium  is  isotropic,  the  tensor  K  becomes  a  scalar  and  G  vanishes. 
In  other  words,  an  irrotational  electric  field  in  an  isotropic  medium  cannot 
induce  a  magnetic  field.  In  an  anisotropic  plasma  G  may  be  expressed  as 
follows : 


G  =  jwe  V  X  K  V  ^ 


=  jW  £  V  X  (K  -  IK  )  V 
o  o 


where  I  is  the  unit,  diagonal  tensor 


This  can  be  expanded  and  simplified  with  the  help  of  the  equation 


V-D  =  V-R  E  =  -V-R  V  =  0 


or 


(Jl 

Ux2 


~2  K  a2 
o  o  o 

+  T~2  +  W  T~2 

o  y  u  z 


) 


4«  =  o 


Thus  the  expression  for  G  becomes 


where 


6  ■  c„  *  sf v 


=  d  - 
o  '  31 


x  = 


j(K’-K  ) 

o 


-j(K’-K  ) 

o 


0 

0 


K’J 


The  elements  of  the  tensor  JL  give  an  indication  of  the  magnitude  of  the 
induced  magnetic  field  H.  For  the  lossless  case: 


K”  = 


XY 


1  -  Y 
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XY 


1  -  Y 


(1-Y2)(1-X) 
K  * 


1  -  (X-Y  > 


It  is  evident  that  the  assumption  of  negligible  magnetic  field  would  be 

2  2 

invalid  in  the  vicinity  of  Y  -  1  (gyroresonance)  and  X-Y  rl  lor  the  lossless 
case.  However,  if  sufficient  damping  were  present  the  terms  of  K  would  be 
bounded.  Thus  the  quasi-static  approximation  is  still  valid  provided  that  the 
frequency  is  low  enough  or  that  the  antenna  is  small  enough  in  terms  of  wave¬ 
lengths  . 

The  presence  of  an  induced  magnetic  field  not  only  restricts  the  appli¬ 
cability  of  the  quasi-static  theory  but  also  suggests  that  unusual  electro¬ 
magnetic  effects  may  occur  in  anisotropic  media.  In  particular,  a  very  small 
antenna  could  become  an  effective  radiator  if  the  induced  magnetic  field  had 
a  phase  and  direction  such  that  outward  power  flow  could  take  place  In  Poyntine's 
theorem,  the  outward  power  flow  through  the  surface  s  is  proportional  to  P 
where 


s 


Power  flow  can  take  place  when  H  has  a  component  pe;  rpendi  cular  to  E  An 
examination  of  ?<  shows  that  E  induces  a  perpendicular  H  through  the  off- 
diagonal  terms  j(Kv-Kq)  Thus  the  magnetic  field  is  90°  out  of  phase  with 
the  electric  field  as  long  as  the  electric  field  has  no  spatial  phase  variation 
(since  H  is  obtained  from  a  volume  integration  of  E  over  all  space).  However, 
if  there  is  spatial  phase  variation  in  the  electric  lield,  the  magnetic  field 
can  have5  an  in-phase  component  so  that  outward  power  flow  is  possible  .  In  2.3.3 
it  will  be  shown  that  an  electric  field  with  spatially  varying  phase  is  produced 
by  a  short  antenna  when  the  differential  equation  1  or  the  field  is  hyperbolic  in 
the  space  co-ordinates  (cf  the  field  formulas  on  p  ige  70) 

The  expression  for  power  flow  P  can  be  put  into  quasi-static  form  by  sub¬ 
stituting  E  -  -  V  +  and  y  x.  H  jeaD  Thus 


65 


P  =  -ju 


n  ds 


The  above  expression  for  P  is  quite  useful;  for  instance  its  analogous  form 

for  the  magneto-static  problem  gives  a  result  identical  to  that  obtained  by 
16 

Trivelpiece  et  al  in  their  Equation  (A, 10). 

2.3,3  The  Differential  Equation  for  the  Field  of  a  Short  Dipole 
If  it  is  assumed  that  H  =  0,  then  E  may  be  expressed  in  terms  of  a  scalar 
potential . 

E  =  -  7  +  (2  .25) 

If  q  is  the  charge  density  in  the  medium  then  V  •  D  =  q  where  D  .=  e  E  =  f^KE. 
This  can  be  written  as 

V  *  K  V  +  =  -  —  (2.26) 

o 


which  is  a  general  expression  of  the  required  differential  equation.  Note 
that  it  reduces  to  Poisson's  equation  when  K  becomes  a  scalar.  In  the  expansion 
of  the  differential  equation  the  off-diagonal  terms  of  K  cancel  (due  to  skew 
symmetry)  and  the  equation  becomes  a  modified  Poisson's  equation  of  the  form 


+  + 

XX  yy 


-b  q 


(2 .27) 


where  1/a  =  K  /K',  b  -  1,  f  K'  and  the  subscripts  indicate  partial  differen- 

o  o 

tiation . 

The  value  of  a  =  \/K'/ K  is  important  since  it  is  a  measure  of  the  degree  of 

o 

anisotropy  of  the  medium.  It  is  necessary  to  define  a  principal  value  of  "a"  to 
facilitate  the  inverse  Fourier  transform  to  be  taken  later;  the  principal  value 
is  defined  such  that  the  argument  of  "a”  is  between  -77/2  and  ■>•77/2.  Since  any 
plasma  medium  is  lossy,  this  guarantees  that  "a"  always  has  a  positive  real  part. 
If  it  is  desired  to  consider  a  "lossless"  plasma  it  is  necessary  in  general  to 
postulate  a  small  loss  and  let  it  approach  zero  in  order  to  select  the  proper 
value  ot  "a" . 

Although  a  lossless  medium  cannot  be  realized  in  practice,  the  concept  of 
such  a  medium  is  definitely  an  aid  to  understanding.  For  a  lossless  anisotropic 
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plasma  K  and  K'  are  both  real  but  may  be  positive  or  negative.  In  contrast 
o 

for  a  uniaxial  crystal  the  values  of  K  and  K’  would  always  be  positive  and 

2  °  2 
therefore  a  would  always  be  positive.  In  the  gyrotropic  plasma  when  a 

is  negative  the  modified  Poisson's  equation  becomes  hyperbolic  whereas 
o 

normally  (a  positive)  it  is  elliptic.  This  has  been  observed  in  the  case  of 
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ferrite  media  by  Damon  and  Eshbach  who  considered  only  the  source-free 

differential  equation.  However  the  real  significance  of  the  elliptic-hyperbolic 

classification  becomes  apparent  only  when  a  source  term,  (-bq),  is  included. 

The  source  is  a  type  of  discontinuity  and  discontinuities  "propagate"  (i„e, 

cannot  vanish)  on  the  characteristic  surfaces  of  a  differential  equation.  The 

characteristic  surface  of  an  elliptic  equation  is  complex  and  thus  has  no 

physical  significance.  On  the  other  hand  the  characteristic  surface  of  a 

hyperbolic  equation  is  real  and  of  considerable  significance. 

Consider  a  point  charge  in  a  "hyperbolic"  medium  (for  the  properties  of 

the  medium  determine  the  differential  equation  classification).  A  charge 

which  goes  to  infinity  at  a  point  will  produce  a  field  which  goes  to  infinity 

everywhere  on  the  characteristic  surface  passing  through  that  point.  Thus 

the  conditions  under  which  the  differential  equation  is  hyperbolic  are  of 

major  importance  in  determining  the  nature  of  the  fields.  These  conditions 

are  presented  graphically  in  Figure  2.5.  Note  in  Figure  2.5  that  the  "operating 

point"  on  the  chart  moves  horizontally  with  a  change  in  electron  density  (as 

in  a  decaying  plasma),  vertically  with  a  change  in  magnetic  field,  and  along 

a  straight  line  through  the  origin  with  a  change  in  frequency. 

The  equation  for  the  characteristic  surface  may  be  obtained  easily 
51 

(see  Sneddon  ,  for  instance).  If  the  modified  Poisson's  equation  in  cylin¬ 
drical  coordinates  is 


^PP  +  V  L 2  *zz  =  -bq 

d 


(2.28) 


then  the  equation  for  the  characteristic  surface  may  be  obtained  as  follows 
(where  the  dot  represents  differentiation  with  respect  to  some  parameter): 


.2  1  -2 
z  +•  -5  P  -  0 


(2.29) 


(PROPORTIONAL  TO  MAGNETIC  FIELD) 
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X  (PROPORTIONAL  TO  ELECTRON  DENSITY)  - ► 


Figure  2.5.  The  elliptic  and  hyperbolic  regions 

Note:  0  is  the  angle  (with  respect  to 
the  Z  axis)  of  the  characteristic  cone 
when  the  differential  equation  is  hyperbolic 
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i 


•  + 
'L  ~  - 


dz 


+ 


(2.30) 


z  =  -  j  P  +  const  .  . 

a 


P 


Under  lossless  hyperbolic  conditions  'V  is  imaginary  and  the  characteristic 
surfaces  are  right  circular  cones  with  half-angle  9  given  by 


tan  0  =  —  -  ~ 
z 


(2.31) 


Thus  an  oscillating  point  charge  in  a  hyperbolic  medium  generates  fields  which 
become  infinite  everywhere  on  the  conical  surface  with  apex  at  the  point  charge. 
Thus  the  most  prominent  feature  of  the  field  can  be  determined  without  carrying 
out  a  detailed  analysis. 

It  is  also  worthwhile  to  see  what  restrictions  the  differential  equation 
imposes  on  a  wave-type  solution. 

If  it  is  assumed  that  q  =  0  and  that  v);  has  the  form  e-^’1*,  then  the 
relation 


(2.32) 


is  obtained. 


This  means  that  k  can  be  real  if  K  /K  is  negative.  Thus  a 
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wave  solution  can  exist  but  the  direction  of  the  wave  normal  k  is  fixed  in 
space.  Its  angle  with  respect  to  the  z-axis  is  <p  which  is  given  by 


tan  (p  = 


(2.33) 


Thus  the  relation  between  0  and  (p  is  given  by 

tan  (p  =  -  cot  0  (2.34) 

where  <p  gives  the  direction  of  phase  progression  and  0  gives  the  direction  in 
which  discontinuities  propagate. 

The  detailed  analysis  of  a  short  dipole  is  carried  out  in  the  following- 
way.  First  a  line  current  distribution  J  is  assumed  and  using  the  equation 
of  continuity 


V-J  +  ju»  q  =  0  ,  (2.35) 

an  equivalent  charge  distribution  q  is  obtained.  Second,  the  potential  4* 
is  determined  using  the  Fourier  transform  pair 


'Y(k)  = 


f(?)  e"jk'r  dr 


(2.36) 


f(r) 


1 

(2  m3 


/v  _ 

f  (k) 


e  dk 


(2.37) 


where  r  =  (x,  y,  z)  and  k  =  (k  ,  k  ,  k  ). 

’  x  y  z 

Third,  the  input  impedance  Z .  ^  is  calculated  approximately  by  integrating 
E-J  over  the  antenna  surface.  The  field  E  is  taken  to  be  the  field  of  a  line 
source  at  the  center  of  the  thin,  cylindrical  dipole,  and  the  current  J  is 
assumed  to  be  spread  over  the  surface  of  the  dipole.  This  method  of  impedance 
calculation  is  essentially  t  lie  same  as  the  "induced  EXIF  method  mentioned 
often  in  the  literature 
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2.3.4  The  Field  and  Impedance  Formulas 

For  the  sake  of  brevity,  only  the  results  of  the  analysis  will  be  pre¬ 
sented.  Refer  to  Figure  2.6  for  a  sketch  of  the  coordinates.  First,  the 
expressions  for  the  electric  field  parallel  to  the  dipole  will  be  given. 

The  magnetic  field  is  always  parallel  to  the  z  axis. 

Dipole  Parallel  to  Magnetic  Field 


„  a  f  2  2  2,"1/2  2  2  2,'1/2 

Ez  =  JTvTe  it"  L  t  (P0  +  a  )  +  (Po  +  a  z2  ) 


-  2(p  2  +  a2z  2) 

o  o 


-1/2 


(2.38) 


where 


z 


2 

1 


(z  +  L)2 


Dipole  length  =  2L 


Dipole  Perpendicular  to  Magnetic  Field 


-1/2  -1/2 

Ex  "  JoT4ff€  k'  L  t  «>l  +  a  *o  >  +  <^2  +  »V> 


2  (p  2  -t-  a2z  2) 
o  o 


-1/2 


(2.39) 


where 


p/  -  (xo  -  L>2  ♦  ,o2 
p/  -  <XQ  .  L)2  x  y‘ 


Note  that  E  and  E  are  taken  at  (x  ,  y  ,  z  )  which  is  the  point  P  in  Figure  2  6 
In  the  following  formulas  for  input  impedance,  it  is  assumed  that  the 
dipole  length  (2L)  is  much  greater  than  its  diameter  (2 p) .  When  the  medium 
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Figure  2.6.  The  co-ordinate  system  lor  the  thin,  cylindrical  dipole 
Note:  dipole  length  =  2L,  dipole  diameter  =  2 p 
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is  isotropic,  these  results  are  in  agreement  with  the  formulas  derived  in  the 
appendix. 

Dipole  Parallel  to  Magnetic  Field 


Z.  = 
in 


jw  7T 


itlnr 


1  t  In  a) 


(2.40) 


n 

where  L  »  -p— r 
I  a  | 

Dipole  Perpendicular  to  Magnetic  Field 

zin  ■  jPirh;-  i  - 1  -  m  i|i]  (2.4i 

where  L  »  | a |  p . 

An  examination  of  Equations  (2.40)  and  (2.41)  shows  that  and  K' 
play  the  parts  of  scale  factors.  Thus  the  two  formulas  can  be  written  as 
one  with  direction-dependent  scale  factors.  The  impedance  formula  with 
normalized  frequency,  length  and  radius  is 


Z  . 
in 


- - — -~c -  [In  —  -  1  ] 

jw  7T  e  L  1  p 
n  on  n 


(2.42) 


where  «  =  «  K'Jk  ,  L  =  hAfK1  ,  p  =  P/J*n. 

n  "  o>  n  L’  n  p 

dependent.  When  the  dipole  is  parallel  to 
K  =  K' .  When  the  dipole  is  perpendicular 
and  K  is  given  by 

r 


Here,  K  and  K  are  direction 
L  p 

the  magnetic  field,  K,  =  K  and 
’  L  o 

to  the  magnetic  field,  K  =  K' 

Lj 


Equation  (2.42)  has  exactly  the  same  form  as  Equations  (2.40)  and  (2.41) 

when  the  medium  is  free  space  Thus  the  free  space  impedance  formula  may 

be  used  if  the  antenna  dimensions  and  the  frequency  are  suitably  scaled.  In 

a  similar  way  electric  field  Equations  (2  38)  and  (2  39)  may  be  scaled  to 

have  the  form  of  fields  in  free  space.  In  the  scaled  formulas  the  electric 

field  components  are  "k/k  E  and  -J  K'  E 

v  o  z  »  x 
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2,3,5  Discussion 

2 

In  the  hyperbolic  region  for  the  lossless  case  (a  is  negative  and  real), 

the  E  field  exhibits  the  conical  discontinuities  predicted  by  the  study  of 

the  characteristic  surfaces  of  the  differential  equation.  When  derived  for 

an  infinitesimal  dipole  the  field  expressions  are  identical  to  the  near 
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field  formulas  obtained  by  Mittra  .  The  expressions  are  similar  to  those 
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of  Kononov,  Rukhadze  and  Solodukhov  ;  however  their  work  predicts  zero 
field  outside  the  characteristic  cone,  a  result  obtained  if  one  does  not 
postulate  a  slight  loss. 

The  expressions  for  impedance  each  have  a  logarithmic  term  containing 
the  anisotropy  factor  "a".  In  a  lossless  medium  under  hyperbolic  conditions 
"a"  is  imaginary  and  therefore  the  logarithmic  terms  will  be  complex.  This 
means  that  the  impedance  has  a  real  part  and  therefore  energy  is  transmitted 
irreversibly  into  the  anisotropic  medium.  At  first  glance  this  result  seems 
very  strange  especially  since  only  quasi-sta i ic,  near-l'ield  theory  was  used. 

However  the  hyperbolic  regions  are  identical  to  the  quasi -static  propagation 
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regions  observed  by  Trivelpiece  and  Gould.  Furthermore,  in  the  hyperbolic 
regions  a  magnetic  field  is  induced  and  electromagnetic  power  flow  can 
take  place  (as  shown  earlier  in  this  chapter).  Tn  the  following  section 
impedance  calculations  are  displayed  graphically  and  compared  with  experi¬ 
mental  results. 

A  major  source  of  theoretical  inaccuracy  may  arise  from  the  method  of 
impedance  calculation  in  which  the  fields  at  the  antenna  surface  are  assumed 
to  be  those  of  a  line  source  of  current  In  an  elliptic  medium  the  fields 
are  continuous  and  this  type  of  approximate  calculation  is  known  to  give  good 
results.  However,  the  discontinuous  fields  in  a  hyperbolic  medium  are  very 
much  dependent  on  the  nature  of  the  assumed  source,  and  thus  the  approximate 
calculation  given  may  be  quite  inaccurate.  A  more  accurate  calculation  might 
result  if  the  antenna  were  considered  to  be  inside  a  small  close-fitting 
free  space  cavity  in  the  anisotropic  medium  (it  is  hoped  that  this  calculation 
will  be  included  in  a  later  report) 

2.4  The  Laboratory  Experiment 

2.4  1  The  Early  Experimental  Apparatus 

The  results  from  this  apparatus  (sec  Figure  27)  were  entirely  qualitative 
but  still  quite  useful  Experimentation  indicated  the  importance  of  the  ion 


Figure  2.7.  Early  experimental  apparatus 
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sheath  and  demonstrated  the  effects  of  RF  heating  of  the  plasma  Experience 
gained  with  the  discharge  tube  and  the  vacuum  system  aided  in  the  design  of 
improved  apparatus  suitable  for  quantitative  measurements. 

DC  bias  was  applied  to  the  probe  using  the  discharge  cathodes  as  refer¬ 
ence  electrodes.  Negative  and  positive  bias  voltages  produced  the  effects 
on  the  reflected  wave  shown  in  Figure  28,  with  apparent  sheath  collapse  at 
+9  volts.  It  is  evident  that  the  impedance  of  the  antenna  is  very  much 
dependent  on  the  size  of  the  ion  sheath  Note  that  Figure  2  8  shows  the 
output  from  a  square-law  detector  as  seen  on  the  oscilloscope  and  thus  the 
vertical  scale  is  proportional  to  RF  power,  not  voltage 

The  effect  of  RF  heating  is  illustrated  in  Figure  29  No  heating  was 
detected  below  1  mm  but  an  appreciable  effect  was  observed  at  20  mw,  the 
maximum  power  available.  The  heating  power  threshold  of  .1  mw  corresponds 
to  a  voltage  threshold  of  .07  volt  (with  a  50  ohm  line).  It  was  expected 
that  the  threshold  would  be  about  equal  to  the  average  energy  of  an  electron 
at  room  temperature,  about  .04  volt  Figure  2  9  shows  also  that  negative 
bias  reduces  the  effect  of  heating  considerably  in  comparison  with  collapsed 
sheath  conditions.  In  addition,  heating  always  produces  the  same  type  of 
change  as  negative  bias  This  suggests  that  heating  produces  sheath  expansion 
as  predicted  by  the  theory  of  Section  2.2.2 

Quenching  ol  the  light  output  of  a  plasma  takes  place  when  the  electrons 
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are  selectively  heated  with  an  RF  voltage  A  photomultiplier  was  located 

in  two  positions  as  shown  in  Figure  2  7  and  power  input  was  20  mw  at  4.3  Gc. 
In  position  A  quenching  of  the  afterglow  was  detected  easily  but  in  position 
B,  quenching  could  not  be  detected  at  all.  This  indicated  that  the  RF  heating 
was  confined  to  the  immediate  vicinity  ol  the  probe. 

Second  harmonic  output  was  mca.su ted  with  20  mw  input  at  4.3  Gc  using 
a  superheterodyne  detection  system  connected  to  the  directional  coupler 
Harmonic  output  was  coniined  t  >  tin  first  lot)  p.s  of  the  afterglow  and  reached 
a  peak  about  20  |is  after  t  hi-  discharge  pulse  Tile  peak  second  harmonic  output 
was  35  db  below  the  fundamental 

The  pulsed  DC  discharge  was  not  •ntirelv  satisfactory  since  the  discharge 
density  often  varied  from  pulse  to  pulsi  as  did  the  balance  between  the  two 
ends  of  the  discharge  tube  This  was  dtu  either  to  momentary  arcs  forming 


Figure  2.8.  Reflected  wave  magnitude  as  a  function  of  probe 
bias.  The  probe  is  matched  with  no  plasma 
(1  <  tQ) .  t  -  tQ  indicates  beginning  of  discharge 


0  t0  TIME ( t )  - ►  t  =  ims. 

Figure  2.9.  The  effect  ol  heating  as  a  function  of  probe  bias 


- RF  input  power  =  . 1  mw 

- RF  input  power  =  20  mw 
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at  irregularities  on  the  discharge  cathodes  or  to  instability  in  the  local 
discharges  (as  evidenced  by  random  shifts  in  the  glow  structure)  Proper 
operation  could  be  obtained  only  by  careful  adjustment  of  the  local  dis¬ 
charge  voltages  and  the  balancing  resistors  (shown  as  100  ohms  in  Figure  2.7) 

As  a  result  of  these  observations  the  following  recommendations  were  made: 

1.  The  discharge  tube  should  be  single-ended,  that  is  it  should  have 
a  single  discharge  cathode. 

2.  The  cylindrical  shape  of  the  discharge  cathode  (15-E  triode  with 
grid  removed)  is  satisfactory,  but  the  cylinder  should  have  a  much 
greater  volume  to  ensure  dependable  starting  of  the  pulsed  discharge 

3.  The  cathode  surface  should  be  smooth  and  its  shape  regular  to  dis¬ 
courage  arcing. 

4.  A  metallic  "ground  plane"  surrounding  the  probe  should  be  provided 

so  that  a  state  of  sheath  collapse  could  be  attained  everywhere 

in  the  vicinity  of  the  probe.  For  good  bias  control  of  the  ion 

sheath,  the  following  area  relationship  is  necessary  (area  of 
cathode)  »  tv'ea  of  "ground  plane")  »  (area  of  probe). 

These  ideas  were  carried  out  successfully  in  the  improved  apparatus  described 
in  the  next  section. 

2.4  2  The  Improved  Vacuum  System  and  Discharge  Tube 

The  vacuum  system  used  for  the  quantitative  experiments  is  shown,  in  Fig¬ 
ure  2.10.  It  proved  to  be  quite  satisfactory  and  a  gas  fill  could  be  left  in 
the  system  for  about  two  days  with  no  apparent  effects  of  leak  contamination. 

_  g 

The  pump-down  period  for  a  pressure  of  10  mm  varied  from  one  to  four  days 
depending  on  such  factors  as  the  amount  of  new  vacuum  grease  in  the  system. 

In  general  it  was  found  advisable  to  ou t gas  the  glass  discharge  tube  with  a 
spark  coil.  The  mam  disadvantages  of  'he  system  were  excessive  grease  (stop 
cocks,  manometer,  0  -ring)  and  slow  pumping  due  to  the  long,  narrow  high- 
vacuum  line. 

The  improved  design  of  the  discharge  cathode  (Figures  2  11  and  2  12) 
resulted  in  a  very  steady,  intense  dischaigc  with  no  arcing  a>  all.  The 
discharge  had  a  tendency  to  concent i ate  at  the  tip  of  the  probe  but  this 
effect  was  greatly  reduced  when  l  lie  magnetic  field  was  applied.  As  was 
mentioned  in  Section  2.2.4,  an  insulating  deposit  1 orms  sm  the  probe  and 
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Figure  2.10  The  vacuum  system 
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must  be  cleaned  off  from  time  to  time.  It  can  be  seen  that  the  probe  dimensions 
and  the  discharge  tube  diameter  are  small  compared  to  a  wavelength.  Therefore, 
this  configuration  is  satisfactory  only  for  testing  the  quasi-static  theory. 

2.4.3  The  Method  of  Impedance  Measurement 

The  impedance  measurements  are  made  using  a  slotted  line  as  shown  in 
Figure  2.11;  a  photograph  of  the  oscillograph  trace  (probe  voltage)  is  taken 
at  each  of  four  points  along  the  slotted  line,  the  points  being  spaced  1/8 
wavelength.  Thus  for  any  desired  time  in  the  afterglow,  a  point  on  the 
Smith  Chart  may  be  obtained  by  purely  graphical  steps  as  follows. 

1.  The  oscilloscope  trace  is  measured  with  dividers. 

2.  This  measurement  is  applied  to  a  parabolic  calibration  chart  which 
compensates  for  the  crystal  detector  square-law  response  and  normalizes 
tne  resulting  voltage  for  the  desired  size  of  Smith  Chart. 

3.  The  resulting  measurement  sets  the  radius  of  an  arc  which  is  drawn 
on  the  Smith  Chart.  Only  three  such  arcs  are  really  necessary  but 
the  extra  one  is  helpful  as  a  check.  Thus  a  point  on  the  chart  can 
be  plotted  quite  accurately. 

In  general  one  must  try  to  minimize  reflections  from  either  end  of  the  slotted 
line,  since  even  small  reflections  can  introduce  considerable  error. 

An  independent  means  of  electron  density  determination  is  not  provided 
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in  these  experiments.  It  is  intended  to  use  the  resonance  probe  method  ’ 
for  this  purpose  in  later  research.  It  is  also  intended  to  provide  a  Langmuir 
probe  system  to  measure  both  electron  temperature  and  density. 

2.4.4  The  Impedance  Results 

The  impedance  results  are  plotted  on  Smith  Charts  (Figures  2.13  to  2.16) 
by  the  graphical  method  described.  The  Smith  Charts  are  normalized  to  50  ohms, 
the  characteristic  impedance  of  the  coaxial  line.  The  collision  frequency  is 
taken  from  Figures  2.3a  and  2.3b  note  that  neon  at  7  mm  and  helium  at  1.2  mm 
have  the  same  electron-molecule  collision  frequency.  Z  varies  from  about  .03 
at  low  electron  densities  to  about  08  at  high  densities  It  was  assumed  in 
making  the  calculations  that  the  plasma  is  m  thermal  equilibrium  at  300°  K. 

The  agreement  between  theory  and  experiments  is  not  very  good  at  low 
magnetic  fields  and  quite  good  at  high  magnetic  fields  By  reducing  transverse 
diffusion,  the  magnetic  field  probably  renders  the  plasma  more  uniform.  For 
zero  magnetic  field  the  experimental  curves  appear  to  suggest  high  energy  loss. 
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Figure  2.13,  Monopole  impedance  for  the  case  Y' 
as  a  Junction  of  electron  density 
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This  could  be  a  result  of  plasma  non- uni formi ty  which  might  increase 
the  apparent  resistance  of  the  plasma  or  permit  some  form  of  surface  wave  to  carry 
energy  away  from  the  probe.  Another  possibility  is  that  the  electron  tempera¬ 
ture  decay  may  have  a  longer  time  constant  than  that  given  by  Dougal  and 

46 

Goldstein  (see  Section  2.2.4,  Equations  (2.23)  and  (2.24))  The  theoretical 

curves  exhibit  a  "kink"  in  the  vicinity  of  plasma  resonance  but  in  the 

experimental  curves  this  kink  is  almost  indistinguishable.  This  smoothing 

effect  could  be  due  to  non-uniformity  in  the  electron  density  or  to  electron 

o 

temperatures  considerably  higher  than  300  K. 

Neon  and  helium  produce  similar  impedance  results  but,  on  the  whole,  the 
neon  curves  are  closer  to  the  theoretical  curves.  This  is  to  be  expected 
since  neon  has  the  higher  recombination  coefficient  and  therefore  should  give 
a  more  uniform  plasma.  Also,  since  the  neon  is  at  a  higher  pressure  than 
the  helium,  the  neon  should  be  less  subject  to  contamination. 

To  summarize,  there  is  good  qualitative  agreement  between  theory  and 
experiment.  It  is  especially  interesting  that  the  impedance  plot  shifts 
from  one  side  of  the  Smith  Chart  to  the  other  as  the  magnetic  field  passes 
through  cyclotron  resonance.  Plasma  non- uni formi t y  and  elevated  temperature 
are  suggested  as  reasons  for  the  poor  quantitative  agreement  between  theory 
and  experiment  especially  evident  at  low  values  of  magnetic  field  The 
problem  of  non-uniformity  is  considered  in  the  following  chapter. 

2.5  The  Effect  of  Non-Uniform  Electron  Density 

In  any  decaying  plasma,  recombination,  ambipolar  diffusion  and  ion  sheath 
formation  will  all  be  present.  Far  from  the  container  walls  volume  recombina¬ 
tion  should  tend  to  dominate  and  maintain  a  uniform  plasma  density.  Closer 
to  the  wall,  ambipolar  diffusion  dominates,  causing  the  electron  and  ion 
densities  to  decrease  Very  closj  to  the  wall  there  is  pr >duced  an  ion  sheath 
in  which  the  electron  density  drops  exponentially.  If  the  wall  is  a  conductor, 
it  may  be  possible  to  collapse  the  sheath  using  positive  bias.  From  these 
considerations  it  is  clear  that  the  plasma  may  be  quite  non-uniform,  especially 
in  the  vicinity  of  a  boundary  such  as  an  antenna  surface 

The  problem  of  non-uniformity  was  considered  by  Tonks^  who  noted  that 
a  paral 1  el -pi  ate  capacitor'  containing  a  non-uniform  plasma  would  appear 
to  have  an  unexpectedly  high  loss  il  the  observer  had  assumed  uniformity 
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Tonks  shows  how  a  broad  plasma  resonance  could  result  from  either  a  non-uniform 
plasma  or  a  lossy,  uniform  plasma  Let  us  now  consider  the  effect  of  collisions 
which  were  not  included  in  the  calculations  reported  by  Tonks  We  shall 
determine  the  impedance  of  an  lsotrop.c  plasma  between  two  parallel  plates 
for  four  electron  density  profiles  as  shown  in  Figure  2  17 

For  the  case  of  unit  spacing  between  plates  of  unit  area,  we  have 


t 


For  example,  consider  the  linear  distribution  Xlx)  X^x  as  shown  in  D  of 
Figure  2.17. 
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The  calculations  using  this  type  :>1  lormula  are  shown  in  Figure  2  17.  It  is 
assumed  that  Z  -  r/w  .  1  and  the  impedance  is  normal  l/.ed  to  give  a  free-space 
impedance  of  -j5  in  ordei  to  allow  comparison  wiih  the  experiments 

It  is  evident  from  an  lnsptction  of  Figure  2  17  that  a  non-uniform  plasma 
with  a  low  collision  frequency  may  bt  l  nd  i  s  t  r  tigu  l  shahl  e  1  r  nn  a  uniform  plasma 
with  a  mui.h  higher  collision  frequency  Fn  rnstamr,  an  investigator  who 
measured  curve  B  would  have  concluded  (assuming  uniform  density)  that  Z  .2 

when  the  correct  value  w.i~,  Z  1  Thus  n;  m- mi  l  f  n  mi  t  v  tends  >o  magnify 
the  apparent  collision  frequency 
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The  impedance  of  a  non-uni iorm,  isotropic 
plasma  between  parallel  plates  as  a  function 
of  peak  electron  density  X()  (0  X()  <  8) 

Note:  0  indicates  X  =  1  and  X  indicates  Xt, 
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If  curve  D  is  compared  with  curve  A,  it  will  be  observed  that  the  point 

of  zero  reactance  is  a  Rood  indicator  of  the  average  electron  density.  To 

put  it  another  way,  the  point  X  -  2  of  curve  D  is  very  close  to  the  resistance 

axis  and  for  this  distribution  the  average  normalized  electron  density  is  X  /2. 

o 

This  means  that  the  impedance  of  a  short  probe  can  give  a  good  estimate  of 
average  electron  density,  but  such  measurements  should  not  be  used  to  estimate 
collision  frequency 

The  curves  of  Figure  2.17  (especially  curve  C)  are  strikingly  similar  to 
the  experimental  curves  of  Figure  2  13  (zero  magnetic  field)  This  suggests 
that  non-uniformity  in  the  plasma  is  the  major  reason  for  the  differences 
between  experiment  and  theory.  The  magnetic  field  reduces  transverse  diffusion, 
so  the  theoretical  and  experimental  results  are  much  closer  in  Figures  2.14, 

2.15  and  2.16.  Nevertheless,  diffusion  (especially  longitudinal)  is  not 
entirely  suppressed  and  the  "kinks"  al  X  !  1  are  not  evident  in  the  experiment 
resul ts 

This  exceedingly  elementary  study  ol'  non-uniformity  suggests  that  losses 
will  be  quite  high  as  long  as  some  portion  of  an  isotropic  plasma  is  in  plasma 
resonance.  Thus  phenomena  such  as  the  ion  sheath  cannot  be  represented  by  a 
''step'1  variation  in  electron  density  near  an  RF  probe.  Instead,  a  continuous 
variation  in  electron  density  should  be  considered  especially  when  the  RF 
signal  is  at  a  frequency  lower  than  the  highest  value  of  plasma  frequency. 

2.6  Planning  for _ the  Gyro- Interact  ion  Experiment 

26  1  Descr ipt l on  of  the  Experiment 

The  ionospheric  gyro-interaction  experiment  is  designed  to  investigate 
the  effects  of  raising  the  electron  temperature  in  the  D'  and  E'  layers  by 
means  of  a  gy ro- frequency  transmitter  carried  by  a  rocket  A  second  signal 
sent  1 rom  the  ground  traverses  the  heated  region  and  is  picked  up  by  a  receiver 
in  the  rocket  so  that  the  received  signal  amplitude  is  a  measure  of  the  effect 
of  heating  Thus  there  are  two  rocket  antennas,  the  receiving  loop  antenna 
(sensing  antenna)  which  is  entirely  inside  the  fibetglass  nose  cone  and  the 
heating  dipole  antenna  which  extends  uitwutd  front  the  sides  of  the  rocket. 

In  addition  there  is  a  small  dipole  antenna  mounted  insult  the  tip  of  the 
nose  cone  this  antenna  is  used  to  monitor  continuously  the  field  of  the 
heating  antenna 
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The  quasi-static  theory  should  be  useful  for  predicting  the  heating  antenna 
impedance  and  also  the  field  strength  at  the  nose  cone  monitoring  antenna.  Both 
of  these  are  measured  m  the  rocket  experiment.  in  the  experiment  the  earth’s 
magnetic  field  acts  in  a  direction  approximately  perpendicular  to  the  heating 
antenna.  Thus  the  theory  for  a  dipole  along  the  x-axis  (magnetic  field  along 
the  z-axis)  should  be  suitable  for  approximate  calculations 

As  described  in  Chapter  1  of  this  report  the  rocket  experiment  has  been 
carried  out.  Unfortunately  the  malfunction  of  the  commutator  switch  caused 
the  impedance  experiment  to  be  omitted  entirely  and  made  the  monitor  probe 
voltage  difficult  to  interpret  (since  the  heating  transmitter  operated  only 
at  its  highest  preset  level).  In  addition,  a  bias  circuit  malfunction  in  the 
monitor  probe  amplifier  cut  off  its  output  for  low  levels  so  that  the  leading 
and  trailing  edges  of  the  monitored  heating  pulse  were  not  recorded.  Thus  it 
is  impossible  to  draw  any  conclusions  about  electron  density  or  temperature 
from  the  monitor  probe  output  recording.  Nevertheless,  calculations  for  the 
expected  heating  antenna  field  and  impedance  are  presented  and  could  be 
utilized  in  further  experiments. 

2.6.2  The  Heating  Field  at  the  Nose  Cone  Tip 

The  E  field  parallel  to  the  antenna  as  given  in  Section  2.3.4  is 


For  the  case  of  a  pickup  antenna  on  the  z-axis,  p^  =  and  p 


0  so  that 


a 


Also,  p^  =  3  05  m  and  z  =  1 , 9U  m 

Pj^  =  9  30  and  z ^  =  3  60  m^. 

If  E^  is  proportional  to  some  normalized  field  E,  we  may  define 


2 

In  free  space  E  =  E  (K'  -  1 ,  a  =1)  Thus  the  theoretical  relative  field 
o 

strength  is  |E/E^,  and  values  will  be  obtained  for  the  appropriate  parameters 

X?  Y,  and  Z,  Since  for  small  Z,  1/K'  =  j2Z/X  the  magnitude  of  the  signal 

received  at  the  nose  cone  tip  could  give  a  good  indication  of  the  collision 

frequency  if  the  electron  density  were  known  Thus  the  presence  of  RF  heating 

(with  constant  electron  density)  should  be  readily  detectable. 

Calculations  of  the  field  strength  E/E  I  are  shown  graphically  in  Fig- 
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ure  2,18  (based  on  a  noon  electron  density  profile  given  by  Waynick  and 
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collision  frequencies  given  by  Nicolet  ),  The  effect  of  heating  is  estimated 
by  performing  some  of  the  calculations  using  a  collision  frequency  increased 
by  a  factor  of  ten.  The  calcul  it  ions  suggest  that  the  nose  cone  field  measure¬ 
ment  should  give  a  good  indication  of  the  magnitude  of  RF  heating  at  altitudes 
between  90  km  and  105  km.  In  general,  this  calculation  indicates  that  antenna 
near  field  measurements  may  be  more  useful  than  impedance  measurements  in 
ionospheric  exploration.  In  particular,  the  near  field  of  a  secondary  antenna 
(mounted  at  right  angles  to  the  heating  antenna)  should  be  a  good  indicator 
of  the  average  heating  effect  in  the  immediate  vicinity  d  the  rocket.  (See  Figure  2-20)  . 
2.6.3  The  Impedance  of  the  Heating  Antenna 

The  input  impedance  of  a  short  dipole  perpendicular  to  the  magnetic  field 
(as  given  in  Section  2.3.4)  is 


The  measured  tree-space  impedance  of  the  rocket  heating  antenna  is  approxi¬ 
mately  -j  5000  ohms.  The  computations  based  on  this  value  of  free  space 
impedance  are  shown  graphically  m  Figure  2  19.  These  are  based  on  the  Waynick 

and  Nicolet  models  and  are  normalized  with  respect  to  R  r  1000  ohms  The 

! ) 

point  chosen  tor  matching  network  design  is  a  compromise  it  compensates 
partially  lor  the  electron  density  and  collision  frequency  but  is  chosen  so 
that  the  reflection  coefficient  always  decreases  with  hea'ing  between  85  and 
95  km. 
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Figure  2.18.  Monitor  antenna  voltage  as  a  function  of  altitude  and  heating 
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2.7  Recommendations  for  Future  Research 

Since  the  quasi-static  theory  appears  to  have  some  merit,  it  should  be 
extended.  For  instance,  the  impedance  calculations  might  become  more  accurate 
if  the  antenna  were  taken  to  be  in  a  close-fitting,  free-space  hole  in  the 
plasma.  This  would  have  the  effect  of  reducing  the  magnitude  of  the  field 
discontinuities  under  hyperbolic  conditions  As  another  example,  the  theoretical 
problem  of  an  antenna  inserted  in  the  end  of  narrow  tube  might  be  considered 
in  order  to  duplicate  the  experimental  conditions  more  precisely. 

The  quasi-static  theory  would  also  be  useful  for  the  calculation  of  mutual 
impedances.  Such  calculations  should  be  relatively  simple,  especially  for  the 
case  where  both  dipoles  are  parallel  to  the  magnetic  field.  Mutual  impedance 
measurements  should  give  better  results  than  self-impedance  measurements  since 
the  former  should  be  less  sensitive  to  ion  sheath  variations.  For  this  reason 
it  is  believed  that  a  "mutual  impedance  probe"  would  be  useful  in  ionospheric 
exploration.  In  addition  it  could  be  tested  in  the  laboratory  using  the 
experimental  arrangment  described  in  this  report. 

A  further  application  of  the  quasi-static  theory  would  be  the  re-entry 
sheath  problem.  The  source  would  be  taken  to  be  inside  a  sphere  of  anisotropic 
plasma  surrounded  by  isotropic  plasma  or  free  space.  Quasi-static  theory  would 
be  used  inside  the  sphere  and  full  wave  theory  outside  it.  Using  this  approach 
the  effect  of  a  locally  strong  magnetic  field  could  be  estimated 

It  should  be  possible  to  develop  a  more  general  impedance  theory  by 
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utilizing  recent  research  on  plasma  conductivity  ’  By  this  means  such 

effects  as  Landau  damping  might  be  taken  into  account.  Non-linear  effects 
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such  as  the  resonance  probe  method  ’  also  deserve  further  study,  both 
theoretically  and  experimentally,  and  possibly  with  the  inclusion  of  magnetic 
field  effects. 

Another  theoretical  approach  would  be  to  consider  the  waves  which  can 
exist  on  a  long,  thin  wire  in  an  anisotropic  plasma.  This  could  assist  in 
understanding  the  behavior  of  antennas  with  dimensions  comparable  to  a  wavelength. 

In  summary,  there  are  many  interesting  and  possibly  fruitful  approaches  to 
the  problem  of  antenna  impedance  in  plasma  media.  In  particular,  ionospheric 
and  laboratory  experiments  to  measure  mutual  impedance  are  strongly  recommended. 

A  proposed  ionospheric  experiment  is  sketched  in  Figure  2  20  The  cyclotron 


Figure  2.20.  A  proposed  rocket  experiment 
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heating  and  sensing  design  is  similar  to  the  design  already  tested,  the  sensing 
wave  transmission  being  used  to  determine  the  heating  effect  at  points  distant 
from  the  rocket.  The  difference  is  in  the  addition  of  an  auxiliary  probe 
antenna  and  its  near-field  monitor.  Measurement  of  the  field  and  impedance 
of  the  probe  antenna  should  permit  fairly  accurate  estimation  of  electron 
density  and  temperature,  and  especially  changes  in  these  parameters. 

2.8  Conclusions 

A  quasi-static  theory  for  the  field  and  impedance  of  a  short  dipole  has 
been  worked  out.  Laboratory  impedance  measurements  for  a  probe  in  an  aniso¬ 
tropic  plasma  agree  fairly  well  with  the  theory,  especially  for  the  case  of 
a  moderate  to  high  magnetic  field  The  theory  predicts  the  excitation  of 
radiating  waves  by  the  short  dipole  but  experiments  have  not  yet  positively 
identified  this  effect 

Differences  between  the  theory  and  the  experiments  probably  are  due  to 
non-uniformity  in  the  plasma  electron  density  which  causes  an  apparent 
increase  in  the  plasma  collision  frequency  If  a  short  antenna  is  used  as 
an  impedance  probe,  plasma  non-uniformity  can  seriously  affect  collision  fre¬ 
quency  measurements,  However,  the  average  electron  density  can  be  determined 
much  more  accurately  (for  instance,  by  observing  the  frequency  at  which  the 
reactance  vanishes). 

The  quasi-static  theory  was  used  to  predict  the  impedance  of  the  heating 
antenna  used  in  the  gyro-interaction  experiments  In  addition  the  field  of 
the  heating  antenna  at  the  nose  cone  tip  was  calculated  (a  small  electric  field 
probe  was  located  there)  This  field  magnitude  was  found  to  be  a  good  indicator 
of  local  plasma  heating  Unfortunately,  a  failure  in  the  rocket  instrumentation 
prevented  the  impedance  and  field  measurements  from  being  carried  out. 

The  formation  of  the  ion  sheath  has  been  studied  and  conditions  for  its 
collapse  have  been  obtained  (using  positive  bias)  In  the  laboratory  experi¬ 
ments  probe  bias  was  found  to  have  a  noticeable  effect  on  impedance  and  thus 
the  sheath  was  collapsed  lor  the  quantitative  measurements  described  in  this 
report.  In  the  experiments  the  ion  sheath  introduced  additional  non-uniformity 
and  thus  increased  the  apparent  RF  losses  in  the  plasma 

Recommendations  for  future  research  were  made  with  special  attention  to 
the  use  of  short,  t ocket -mounted  antennas  to  measure  electron  density  and 
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temperature.  Such  measurements  close  to  the  rocket  would  help  in  the  inter¬ 
pretation  of  the  sensing  wave  amplitude,  the  only  recorded  parameter  related 
to  heating  at  some  distance  from  the  rocket.  For  future  gyro-interaction 
experiments,  measurements  of  the  field  and  impedance  of  an  auxiliary  antenna 
perpendicular  to  the  heating  antenna  should  give  a  good  estimate  of  ionospheric 
properties  close  to  the  rocket.  Such  measurements  would  be  especially  useful 
in  studying  changes  in  electron  density  and  temperature  (i.e.  the  transient 
behavior  of  the  ionosphere). 
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APPENDIX  TO  CHAPTER  2 

IMPEDANCE  OF  AN  ANTENNA  IN  A  CONDUCTING  MEDIUM 
Georges  A.  Deschamps 

The  dependence  of  the  input  impedance  of  an  antenna  on  the  properties  of 
the  surrounding  medium  is  particularly  simple  when  the  antenna  is  made  of 
perfect  conductors  and  is  fed  from  a  region  small  in  terms  of  the  operating 
wavelength . 

The  result  can  be  expressed  as  follows;  the  impedance  normalized  to  that 
of  the  medium  takes  the  same  value  at  frequency  w  m  a  medium  of  index  n  as  it 
does  at  frequency  ncj  in  a  medium  of  index  1 

This  is  a  known  result  in  the  theory  of  antenna  modeling  but  it  is  not 
widely  recognized  that  the  relation  holds  when  the  medium  characteristic 
constants  are  complex  numbers,  as  for  example  in  the  case  of  a  lossy  medium. 

The  following  short  proof  will  convince  the  reader  that  this  is  true. 

The  surrounding  medium  is  assumed  homogeneous  and  characterized  by  its  per¬ 
mittivity  and  permeability  constants  £  and  p  where  £  and  p  may  be  complex 

numbers.  Consider  at  frequency  u/  in  the  medium  (£  p  )  the  field  E  (r) 

o  o  ^o  o  ’ 

HQ(r)  which  results  from  applying  one  ampere  current  source  at  the  terminals 

of  the  antenna.  This  field  is  determined  by  Maxwell' s  equations  and  the  proper 

boundary  and  radiation  conditions.  When  the  medium  is  changed  to  (£,  p)  a 

field  which  satisfies  the  new  Maxwell  s  equations  and  the  same  boundary  and 

radiation  conditions  is  obtained  by  scaling  the  vector  E  ,  H  to  E  =  E  (£  /C)1^2 
1/2  o  o  o  o  ’ 

H  =  H^(p^/p)  and  simultaneously  changing  the  operating  frequency  to  w  - 
1  /2 

,•  Integrating  in  the  feed  region  from  one  conductor  to  the 
°  o  o  g 

other  (Figure  A.l)  gives  the  terminal  voltage  V  =  f  E.ds  and  integrating  along 
a  loop  c  around  the  conductors  gives  the  input  current  I  =  J H  ds  Thus 
keeping  these  integration  paths  fixed  when  the  medium  is  changed  we  see  that 
the  impedance  V/I  behaves  in  the  same  manner  as  the  ratio  E/H  at  any  point. 

If  we  denote  by  Z(oj,  £,  p)  the  impedance  of  the  antenna  at  frequency  ^  in  the 
medium  (£ ,  p)  we  have  the  relation 


Z  (uq  £,  p) 


,  >1/2  /  £  \1//2  f  1/2 

w  id  ‘te) . 


(A.l  ) 


Figure  A. 


V 
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Two  conductor  antenna 
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Figure  A. 3.  Sketch  of 


Figure  A. 2.  Antenna  with  TEM  feed 
waveguide 
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Introducing  the  index  n  =r  (u.£/u.  £  )  and  the  intrinsic  impedances  of  the 

1/2  °  °  1/2 

medium  £  =  (u  /£  )  and  £  =  (u./£  )  this  becomes 

o  'o  o  r 

\  Z  (oo,  £,  pj  =  Z(nu),  eQ,  |xo) 
o 

which  can  be  expressed  by  the  statement  made  above. 

If  only  the  permittivity  is  changed -from  €  to  £  = 
simplifies  to 

Z  (oj  €  )  =  -  Z(na),  £  ) 

The  proof  is  obviously  valid  whether  the  various  characteristic  constants 
are  real  or  complex.  The  assumption  about  a  perfect  conducting  antenna  is 
used  to  assert  that,  after  scaling,  the  field  still  satisfies  the  boundary 
condition  (E^  ^  =0.)  f  I f  another  value  of  surface  impedance  was  imposed 
on  the  antenna  it  would  have  to  be  changed,  together  with  the  medium,  so  as  to 
keep  the  same  relative  value  with  respect  to  the  surrounding  medium.]  The 
condition  on  the  size  of  the  source  or  feed  region  is  used  in  identifying 
the  integral  of  the  fields  along  some  paths  with  the  voltage  and  current 
inputs.  Even  with  these  restrictions,  a  large  number  of  useful  antennas 
are  included:  the  dipole  fed  from  a  small  gap,  the  biconical  or  discone 
antenna,  log-periodic  and  log-spiral  antennas  fed  from  a  point 

Strictly  speaking  the  proof  dues  not  apply  to  an  aperture  antenna  because 
one  cannot  find  integration  paths  ab  and  C  in  the  medium  such  that  the  integral 
of  E  and  H  represent  respectively  V  and  I.  The  proof  can  however  be  extended 
to  an  antenna  fed  by  means  oi  a  two  conductor  waveguide  supporting  a  TEM  mode 

provided  the  change  of  medium  from  ('€  ,  p  )  to  (€  o)  is  carried  out  inside 

o  o 

the  waveguide  all  the  way  to  a  terminal  plane  where  the  field  in  regions  I 
and  II  (Figure  A. 2)  can  be  represented  by  the  principal  modes  only. 

An  extreme  example  of  this  is  a  plane  wave,  normally  incident  on  a  slab 
of  material  (£,  p.)  backed  by  a  conducting  plane  If  the  thickness  of  slab 
is  a,  the  impedance  Z  =  E/II  at  the  interface  is 


n  €  the  relation 
o 


(A. 2) 


Z  =  -,j  4  tan  ka 
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an  expression  which  satisfies  Equation  (A.l)  when  the  medium  characteristics 
are  changed.  It  is  easy  to  adapt  the  proof  given  above  to  the  general  TEM 
waveguide . 

When  the  medium  is  modified  only  outside  of  the  actual  aperture  of  the 

antenna,  the  field  in  that  aperture  usually  contains  higher  order  modes  and 

the  relation  (A.l)  may  still  be  applied  but  only  as  an  approximation. 

The  relation  can  be  used  to  find  the  impedance  of  an  antenna  in  a  lossy 

medium,  i.e.,  a  medium  characterized  by  a  complex  value  of  n.  This  requires 

a  knowledge  of  Z(w)  for  complex  values  of  u>.  The  simplest  case  is  when  Z 

is  known  analytically.  Take  for  example  a  short  antenna,  represented  by  a 

capacitance  C  in  free  space  Z  =  1/JCuj.  Changing  the  medium  to  (C ,  p) 

Equation  (A.l)  reduces  to  the  obvious  result,  Z  =  €  /€  1/jCw  If  the  medium 

2  2  2  ° 

is  a  plasma  with  refractive  index  n  =  1  -  to  /to  relation  (A.l),  applied 
to  the  admittance  gives 


Y  (to,  n)  =  n  Y  (nto,  1)  =  jCtu 


I 

jLto 


with 


L  = 


As  a  next  approximation  the  impedance  of  a  short  antenna  takes  the  form  Z  = 
2 

a  to  +  jb/to  and  therefore  in  a  plasma 


„  -if  2  2  ,  2  2 

Z  =  a  to  \/to  -  to  -  lbto  (to  -  w  ) 

VP  P 

When  the  analytic  expression  of  Z  is  not  known  but  Z  has  been  measured 
for  real  values  of  to,  its  values  for  complex  frequencies  (with  negative 
imaginary  part)  can  be  deduced  by  analytic  continuation.  This  can  be  carried 
out  by  graphical  methods  based  on  the  fact  that  the  mapping  w  — »■  Z  must  be 
conformal.  An  illustration  is  shown  in  Figure  A  3  Another  approach  is  to 
note  that  for  negative  values  of  Im  to,  z  is  an  analytic  function  without  singu¬ 
larities,  thus  the  Cauchy  formula  can  be  used 


Z(w)  = 


-r 


where  Z(£)  is  known  on  the  real  axis. 

Since  an  experimental  knowledge  of  Z  is  usually  limited  in  range  and 

accuracy,  there  will  be  serious  limitations  to  the  validity  of  this  procedure 

far  away  from  the  measured  data.  Using  an  analytical  form  which  represents 

the  data  over  a  limited  range  also  permits  continuation  only  in  the  neighborhood 

of  that  range.  However,  for  the  case  of  moderate  losses  (say  for  Im  n  <  10  Re  n) 

even  a  graphical  extrapolation  based  on  a  free  hand  sketch  is  adequate. 

To  test  the  method  and  illustrate  this  last  point  some  results  obtained 

recently  by  R.  W.P.  King  and  C.  W.  Harrison,  Jr.5®were  used.  These  authors 

consider  cylindrical  antennas  of  total  length  2h  and  given  diameter-to-length 

2 

ratio  (1/75).  The  medium  is  characterized  by  p  =  p  6  =  €  -  jCT/co  or  n  = 

e  -  ja/ooe  . 

r  o 

For  these  antennas  let  us  designate  the  input  impedance  at  frequency  u 

by  Z(2h,  to,  n).  Thus  letting  co^  be  the  frequency  and  \  the  wavelength 

corresponding  to  6  me  Table  1  or  Figure  4  of  the  reference  gives  Z(\  /2\if  ,  to  n) 

o  '  r  o' 

for  several  values  of  n.  Thus  only  antennas  whose  length  2h  has  been  adjusted 

to  remain  equal  to  half  of  a  wavelength  in  the  medium  are  considered 

A  change  in  size,  the  medium  being  fixed,  is  equivalent  to  a  change  in 

frequency,  thus  the  expression  is  also  Z(\  /2,  u)  /v/6-,  n)  and  from  relation 

o  o  ’  r 

(A.l)  this  is  1/n  Z(\  /2  n /j€  uj  1)  The  latter  is  the  impedance  of  a 

o  V  r  o 

half-wave  antenna  in  air  for  the  complex  frequency  n^  /\/t  (Note  that  in  all 

o  V  r 

these  formulas  n  must  be  considered  as  fixed  at  the  value  n(u*  )  For  real 

o 

frequencies  we  can  sketch  the  locus  of  Z  from  the  data  in  the  same  paper  We 

can  then  extend  this  locus  to  complex  values  of  w  (with  negative  imaginary 

part)  and  deduce  from  the  plot  values  of  Z(\  /2^~,  ^  .  n)  for  several  CT/uf  . 

o  !  r  o  o 


(from  reference) 


83  +  j  40 
100  *  j  36 
116  +  j  32 
131  4  j  30 
167  t  j  28 


(from  analytic  continuation) 


83  ♦  j  40 


100  -  j  35 
117  i-  j  31 
132*  j  27 
168  *  i  24 


0 
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The  agreement  is  very  good  in  that  range.  For  larger  values  of  0/oj6  it  becomes 

o 

doubtful  and  the  method  of  continuation  would  have  to  be  refined  to  take  into 
account  a  larger  part  of  the  locus  Z  (w)  for  aj  real. 

The  total  complex  plot  of  Z(0j)  suggests  a  simple  method  for  interpreting 
the  result  of  an  impedance  measurement  in  terms  of  the  medium  characteristics. 
The  accuracy  of  this  plot  can  be  improved  by  taking  a  few  measurements  in  a 

lossy  medium.  Making  use  of  the  analyticity  of  ca - >Z  should  help  organize 

the  results  and  reduce  appreciably  the  experimental  effort. 

This  work  was  done  in  connection  with  contract  AF19(604)-5565  The 
computations  reported  here  were  carried  out  by  Mr.  K.  G.  Balmain. 
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3.  PLASMA  NONLINEARITIES  -  J.  T  Verdeyen 


3.1  Introducti on 

The  nonlinear  response  of  plasmas  to  electromagnetic  waves,  as  evidenced 

by  the  generation  of  power  at  the  second  harmonic  of  the  applied  frequency, 

was  studied  during  the  period  of  this  contract.  The  results  of  the  theoretical 

analysis  and  the  laboratory  experimental  study  have  been  given  in  a  separate 
58 

report  .  This  section  will  deal  only  with  the  pertinent  features  of  the 
nonlinear  study  which  apply  to  the  rocket  probe  of  the  ionosphere 

It  was  found  in  the  experimental  study  that  mixing  of  two  electromagnetic 
waves  of  frequencies  f^  and  f  occurred  in  a  plasma  provided  the  amplitude  of 
one  of  the  waves  was  sufficiently  great  It  was  found  that  if  this  amplitude 
was  great  enough  to  cause  an  additional  degree  of  ionization  (i.e.  breakdown) 
then  the  conversion  efficiency  was  greatly  enhanced  over  the  case  without 
breakdown. 

It  was  estimated  that  the  heating  wave  in  the  rocket  would  cause  a  local 
breakdown  in  the  E  layer  of  the  ionosphere,  and  thus  could  perform  the  function 
of  a  local  oscillator  in  a  superheterodyne  receiver.  One  could  then  attempt 
to  detect  one  of  the  products  of  the  mixing  of  the  heating  wave  and  the 
sensing  wave  in  this  localized  region  of  the  ionosphere. 

The  difference  frequency  was  chosen  to  be  the  signal  to  be  detected  in  this 
rocket  probe  for  two  reasons 

1.  The  lower  regions  of  the  ionosphere  would  provide  an  effective  shield 
against  any  spurious  signals  from  other  radio  stations.  Thus,  if 
any  signal  was  detected,  one  could  be  sure  that  it  was  a  product  of 
the  nonlinearity  in  the  ionosphere,  and  not  a  signal  from  the  various 
local  radio  stations. 

2.  The  use  of  the  low  frequency  signal  permits  one  to  use  less  sophisticated 
circuitry  in  the  RF  section  of  the  receiver.  The  weight  and  space 
requirements  of  this  receiver  dictated  that  transistors  be  used 
throughout  The  stability  of  transistor  radio  receiver  is  easier 

to  achieve  at  the  lower  frequencies.  Furthermore  the  difference 
between  the  sensing  and  the  heating  frequency  (610  kc)  was  low 
enough  so  that  a TRF  receiver  was  feasible  rather  than  a  heterodyne 
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scheme.  The  use  of  the  TRF  receiver  avoids  one  source  of  noise  in  the 
heterodyne  scheme — the  first  detector. 

It  should  be  emphasized  here  that  this  part  of  the  rocket  probe  was  con¬ 
ceived  only  four  months  prior  to  actual  firing  of  the  rocket,  and  thus  time 
limitations  and  availability  of  components  dictated  the  choise  of  the  circuitry. 
3.2  System  Concept  of  the  Difference  Frequency  Receiver  (610  kc) 

The  operation  of  difference  frequency  receiver  and  its  place  in  the  rocket 
system  is  best  illustrated  by  the  use  of  the  timing  diagram  in  Figure  3.1. 

The  heating  wave  transmitter  in  the  rocket  was  pulsed  on  for  500  (iseconds  with 

13.3  milliseconds  between  pulses.  The  power  level  of  the  transmitter  was  to 
be  changed  in  ratios  of  1:10:100,  with  each  power  level  being  held  fixed 
during  45  pulses.  This  is  shown  in  Figure  3.1a.  The  sensing  wave  was  broad¬ 
cast  from  the  ground  to  the  rocket  as  a  CW  wave.  It  is  obvious  therefore  that 
a  610  kc  signal  would  only  be  present  during  the  time  interval  of  the  heating 
wave,  and  thus  the  receiver  for  this  610  kc  signal  need  only  be  on  during  this 
time  interval.  The  block  diagram  for  this  receiver  is  shown  in  Figure  3.2. 

The  610  kc  signal  was  to  be  amplified  by  the  TRF  section  and  then  fed  into 
the  gated  clamp.  This  clamp  is  open  for  a  time  duration  (rd  1  msec)  determined 
by  the  monostable  multivibrator,  which  is  triggered  by  the  same  pulse  which 
triggers  the  transmitter.  During  the  duration  when  the  clamp  is  open,  the 
610  kc  passes  through  to  the  detector.  The  output  of  the  detector  is  amplified, 
and  the  energy  contained  in  the  pulse  is  stored  on  a  large  capacitor.  Thus  the 
voltage  on  the  condenser  is  proportional  to  the  energy  contained  in  the  pulses 
of  610  kc  signal.  The  peak  voltage  on  the  capacitor  is  limited  to  5  volts 
by  a  Zener  diode  across  the  output. 

The  voltage  on  this  capacitor  was  to  be  sampled  three  times  during  each 
power  level  of  the  transmitter.  At  the  end  of  each  power  level,  the  capacitor 
was  shorted  to  ground  by  a  reset  switch.  The  signal  for  a  reset  was  derived 
by  a  coincidence  of  two  pulses,  reset  1  and  2,  which  were  fed  into  the  "and" 
gate. 

3.3  Conclusions 

It  has  already  been  noted  that  considerable  data  was  lost  when  the 
telemetering  switch  failed.  Unfortunately,  the  data  from  this  part  was 
among  that  which  was  lost.  However  it  is  fel  t that  any  future  probes  such 
as  this  should  be  instrumented  to  observe  coherait  nonl  mean  ties  such  as 
described  above 


45  PULSES 


Figure  3,1b.  Detail  of  envelope  of  the  heating  wave 
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